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I.  INTRODUCTION 


The  United  States  military  has  undergone  numerous  changes  since  the  end  of  the 
Cold  War.  Specifically,  the  U.S.  Navy  experienced  a  shift  in  the  area  of  engagement 
from  the  “blue  water”  (water  depth  greater  than  100  m)  Soviet  threat,  to  littoral  regions  of 
the  world.  It  is  important  to  predict  instantaneous  currents  and  tides  in  littoral  zones  such 
as  the  continental  shelf,  estuary,  and  rivers.  The  Naval  Oceanographic  Office 
(NAVOCEANO)  implements  and  runs  several  types  of  two-dimensional  (2D)  and  three- 
dimensional  (3D)  circulation  models  for  many  geographical  domains  around  the  world. 
Most  of  these  are  run  in  an  operational  mode  that  generates  daily  products  for  the  Fleet. 
These  products  support  a  wide  variety  of  applications  ranging  from  Special  Operations, 
Mine  Warfare,  Expeditionary  Warfare,  Object  Drift,  Search  and  Rescue,  Chemical  and 
Oil  Spill. 

Complexity  of  coastal  bay  physics  as  well  as  the  forcing  functions  at  the  seaward 
boundary  to  the  open  coast  often  warrants  the  use  of  high  resolution,  three  dimensional 
(3D),  baroclinic  models,  which  are  nested  to  offshore,  data  assimilating  circulation 
models.  The  purpose  of  this  study  is  to  detennine  how  well  2D  depth-integrated  models 
will  satisfy  certain  Navy  applications  in  coastal  bays.  In  regions  dominated  by  tides  and 
local  wind  forcing,  this  assumption  is  reasonable  for  time  scales  ranging  from 
“instantaneous”  (what  a  diver  will  feel  in  the  water)  to  days  (the  drift  of  floating  objects). 
The  effort  to  implement,  tune,  and  validate  these  types  of  models  for  limited,  specific 
applications,  is  a  fraction  of  the  time  than  for  baroclinic,  nested  models. 

This  thesis  provides  guidance  for  identifying  the  baroclinicity  and  detennining 
which  forcing  function  should  be  ignored  or  included  in  modeling  the  littoral  area  of 
interest.  The  San  Diego  Bay  is  selected  for  the  study  because  of  the  plethora  of  data 
available,  its  unique  pure  tidal  forcing  and  its  significance  for  the  U.S.  Navy.  The 
hydrodynamic  model  “Water  Quality  Management  and  Analysis  Package”  (WQMAP)  is 
used  to  predict  the  instantaneous  currents  with  various  forcing  functions  (tides,  winds, 
and  lateral  boundary  fluxes)  and  the  hydrochemical  model  “Chemical  Management  and 

Analysis  Package”  (CHEMMAP)  is  used  to  predict  the  water  contamination  and  simulate 
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chemical  attacks/accidents  in  San  Diego  Bay,  which  raise  considerations  regarding  public 
health,  the  economy,  ecology  or  even  national  security.  Both  models  were  purchased  by 
U.S.  NAVOCEANO.  Finally,  the  study  shows  the  barotropic  nature  of  San  Diego  Bay, 
the  slight  significance  of  wind  and  the  vulnerability  of  a  semi-enclosed  tidal  basin  in  a 
possible  chemical  attack  or  accident. 

The  rest  of  the  thesis  is  organized  as  follows.  Chapter  II  examines  the  geography, 
the  topography  and  the  forcing  in  San  Diego  Bay  in  a  detailed  study  of  the  site.  Water 
exchange  and  water  quality  monitoring  will  also  be  addressed.  Chapter  III  investigates 
the  baroclinicity  of  the  San  Diego  Bay  using  available  data  and  shows  how  and  where 
this  barotropic  nature  cannot  be  assumed  or  applied.  Chapter  IV  describes  the  WQMAP 
model,  its  physics,  possibilities,  assumptions  and  shortcomings.  Chapter  V  evaluates  the 
WQMAP  model.  The  evaluation  will  be  done  through  statistical  analysis  of  u  and  v 
components  of  current  speed  based  on  verified  data  (ADCP  measurements)  but  also 
through  harmonic  analysis  and  decomposition  of  the  tides  produced  by  the  model  and 
comparison  with  verified  data.  Chapter  VI  briefly  describes  the  CHEMMAP  model  and 
presents  the  chemical  elements  used  in  the  study,  explains  the  selections  of  chemicals, 
and  possible  outcomes.  Chapter  VII  presents  the  results  of  the  six  threat  scenarios  and 
their  combinations  to  obtain  general  conclusions  but  also  some  distinct  differences. 
Chapter  VIII  presents  the  overall  conclusions  and  suggests  further  research  and 
applications  for  future  studies. 
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II.  SAN  DIEGO  BAY  OCEANOGRAPHY 


A.  GEOGRAPHY  AND  BOTTOM  TOPOGRAPHY 

San  Diego  Bay  (Figure  1)  is  located  near  the  west  coast  of  southern  California.  It 
is  a  relatively  small  basini  (43-57  km“)  nearly  25  km  long  and  1-4  km  wide.  It  shapes  a 
flipped  “T“  and  extends  to  the  north  to  the  city  of  San  Diego  and  to  the  south  to 
Coronado  Island  and  Silver  Strand,  with  a  northwest  to  southeast  orientation.  It  is  a 
natural  harbor  sheltered  by  overlapping  peninsulas  (in  west  Point  Loma  and  in  east 
Coronado).  San  Diego,  being  the  last  southwest  major  city  neighboring  Mexico,  is  a  very 
important  port,  especially  for  the  U.S.  Navy,  since  it  hosts  the  headquarters  of  the  11th 
U.S.  Naval  District  and  a  large  portion  of  the  American  Fleet. 


Figure  1.  Main  geographical  locations  in  San  Diego  Bay  I  (From: 
http://sdbay.sdsc.edu/html,  last  accessed  on  May  25,  2005). 


1  The  area,  43krrr.  is  measured  at  mean  lower  low  water  (MLLW). 
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|  16.25 ->17.5 
I  17.5 ->18.75 

■  18.75  •>  50 


Figure  2.  San  Diego  Bay  Bathymetry  (From:  WQMAP). 


The  topography  is  not  homogeneous  (Figure  2),  and  the  average  depth  is  of  6.5  m 
(measured  from  the  mean  sea  level).  The  northem/outer  part  of  the  bay  is  narrower  (1-2 
km  wide)  and  deeper  (reaching  depth  of  15  m)  and  the  southern/inner  part  is  wider  (2-4 
km  wide)  and  shallower  (depth  less  than  5  m). 

Since  San  Diego  is  a  semi-closed  bay,  it  exchanges  with  the  Pacific  Ocean  only 
through  a  single  channel  at  the  mouth.  Near  the  mouth  of  the  bay,  the  north-south  channel 
is  about  1.2  km  wide,  bounded  by  Point  Loma  to  the  west  and  Zuniga  jetty  to  the  east 
with  depths  between  7  and  15  m  (Chadwick  -  Largier,  1999a).  The  western  side  of  the 
channel  is  shallower  than  the  east  side. 

B.  TIDES  AND  TIDAL  CURRENTS 

San  Diego  Bay  is  a  perfect  example  of  a  tidal  basin  connected  to  the  ocean  by  an 
inlet  with  an  artificial  jetty  (Zuniga)  built  to  control  beach  erosion.  The  Zuniga  jetty 
extends  almost  one  mile  offshore  Zuniga  Point  and  most  of  it  is  not  clearly  visible  at  high 
water  (Figure  3).  Obviously,  the  bay  has  been  intensively  engineered  to  accommodate 
shipping  activities.  Ninety  percent  of  all  available  marsh  lands  and  fifty  percent  of  all 
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available  inter-tidal  lands  have  been  reclaimed  and  dredging  activities  within  the  bay 
have  been  equally  extensive  (Peeling,  1975;  Wang  et  ah,  1998).  Kelp  forests  extend 
approximately  2  km  south  of  Point  Loma  (Figure  4)  and  along  its  western  side.  They  are 
quite  thick  and  they  create  seasonal  dumping  of  currents  to  about  one-third  their  values 
outside  (Jackson  and  Winant,  1983). 

Currents  in  San  Diego  Bay  are  predominately  produced  by  tides  (Wang  et  al., 
1998).  This  tidal  exchange  between  the  ocean  and  the  bay  is  a  result  of  a  phenomenon 
called  “tidal  pumping”  (Fischer  et  al,  1979).  The  “pumping”  of  water  is  caused  due  to  the 
flow  difference  between  the  ebb  and  the  flood  flows. 

Being  located  at  mid-latitude,  tides  and  currents  within  the  San  Diego  Bay  are 
dominated  by  a  mixed  diurnal-semidiurnal  component2  (Peeling,  1975).  The  tidal  range 
from  mean  lower-level  water  (MLLW)  to  mean  higher-high  water  (MHHW)  is  1.7  m 
with  extreme  tidal  ranges  close  to  3  m  (Chadwick  and  Largier,  1999a).  Typical  tidal 
current  speeds  range  between  0.3-0. 5  m/s  near  the  inlet  and  0. 1-0.2  m/s  in  the  southern 
region  of  the  bay  (Wang  et  al.,  1998).  The  phase  propagation  suggests  that  the  tides 
behave  almost  as  standing  waves  with  typical  lags  between  the  mouth  and  the  back 
portion  of  the  bay  of  10  min  and  a  slight  increase  in  tidal  amplitude  in  the  inner  bay 
compared  to  the  outer  bay.  The  overall  tidal  prism  for  the  bay  is  5.5  x  10  m  and  the  tidal 
excursion  is  larger  than  the  mouth  with  a  value  of  4.4  km  (Chadwick  and  Largier,  1999b). 
Offshore  of  the  bay,  the  average  current  speed  is  0.1  m/s.  The  currents  are  equatorward  in 
all  seasons  and  ninety  percent  of  their  measurements  range  between  1  and  0.25  m/s 
(Jackson  and  Winant,  1983). 

The  form  ratio  (ratio  between  diurnal  -Kl+Ol-  and  semidiurnal  -M2+S2-)  shows 
that  the  tidal  constituents  are  mixed  (Table  1).  Although  there  are  21  harmonic  tidal 
constituents,  it  is  possible  to  approach  the  solution  by  using  four  or  five  (including  N2). 


2  Hence,  time  series  data  collection  should  cover  at  least  two  weeks,  preferably  32  days  to  account  for 
spring-neap  variations  (Wang  et  al.,  1998) 
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Figure  3.  Zuniga  point  and  jetty  in  the  entrance  of  San  Diego  Bay. 


Figure  4.  Point  Loma  in  the  western  entrance  of  San  Diego  Bay. 
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Name 

Amplitude  (m) 

Epoch  (degrees) 

1 

M2 

0.576 

148.9 

2 

S2 

0.233 

145.9 

3 

N2 

0.136 

128.7 

4 

K1 

0.352 

210.5 

5 

01 

0223 

195.6 

6 

NU2 

0.027 

134.3 

7 

MU2 

0.010 

109.7 

8 

2N2 

0.018 

108.7 

9 

OOl 

0.010 

225.4 

10 

LAM2 

0.004 

147.5 

11 

Ml 

0.011 

194.2 

12 

J1 

0.018 

217.9 

13 

SSA 

0.017 

272.7 

14 

SA 

0.063 

182.0 

15 

RHO 

0.008 

189.2 

16 

01 

0.041 

188.7 

17 

T2 

0.014 

145.9 

18 

2Q1 

0.006 

180.7 

19 

PI 

0.109 

208.8 

20 

L2 

0.013 

121.7 

21 

K2 

0.065 

139.3 

Table  1 .  The  mixed  diurnal-semidiurnal  nature  of  San  Diego  Bay  tides.  (From: 

National  Ocean  Service  (NOS)  accepted  harmonic  constants  for  station 
number  9410170  in  San  Diego,  CA  Latitude:  32°  42.8’  N  Longitude: 

117°10.4’W). 


C.  WATER  EXCHANGE 

Since  San  Diego  Bay  has  a  relatively  narrow  mouth  and  a  large  portion  of  shallow 
water,  the  percentage  of  total  bay  water  exchanged  during  a  tidal  cycle  can  be  quite 
significant.  At  mean  low  low  water,  the  bay  has  a  surface  area  of  4.3x10  m  and  a 
volume  of  2.795  x  10  m  .  At  the  mouth  of  the  bay,  the  tidal  prism,  or  volumetric  flux 
passing  a  cross-section  during  a  single  flooding  cycle,  can  reach  40%  of  the  mean  volume 
of  the  bay  during  a  strong  flood.  For  an  average  tidal  cycle,  the  tidal  prism  at  the  mouth  is 
about  13%  of  the  mean  total  bay  volume.  During  a  single  tidal  cycle,  water  from  the  front 
portion  of  the  bay  mixes  with  open  ocean  water  and  exchanges  with  the  bay  water  that 
existed  on  the  previous  cycle.  This  open  ocean  water  exchanges  further  into  the  bay  on 
successive  cycles  (Figure  5). 
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Time  for  exchange  of  50%  of  total  bay  volume  with  the  ocean  water 


Time  (Hours) 

Figure  5.  Time  for  exchange  of  50%  of  total  bay  volume  with  the  ocean  water. 
(From:  http://sdbay.sdsc.edu/html/modeling2.htmh  last  accessed  on  May  25, 

2005). 


D.  HYDROGRAPHIC  CONDITIONS 

Freshwater  inflow  to  the  bay  is  minimal,  comes  less  from  the  Otay  River  and 
mainly  from  the  Sweetwater  River  located  to  the  southern  part  of  the  Bay  and  occurs  only 
during  winter  storms.  Both  rivers  are  regulated  by  storage  reservoirs.  The  San  Diego 
River  has  been  diverted  by  the  U.  S.  Army  Corps  of  Engineers  since  1875  and  no  longer 
empties  into  San  Diego  Bay. 

San  Diego  Bay  can  be  regarded  as  a  vertically  well-mixed  estuary  (Wang  et  ah, 
1998).  In  the  south  bay,  the  currents  are  much  smaller  than  in  the  north  and  the  model 
performance  cannot  be  well  documented.  The  average  water  temperature  in  San  Diego 
Bay  is  21°C  and  ranging  between  14°C  and  26°C.  The  average  temperature  during  the 
summer  (late  June  to  late  August)  is  23°C.  The  temporal  variations  of  temperature  can  be 
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seen  in  Figure  6.  Using  data  from  five  different  stations  of  the  Bay  Wide  Water  Quality 
Monitoring  Program  of  the  Port  of  San  Diego,  the  conclusion  was  that  there  is  zonal 
separation  as  regards  the  water  temperature. 


Figure  6.  Water  temperature  variation  in  San  Diego  Bay  from  8/28/93  to  8/28/94 
(From:  http://www.co-ops.nos.noaa.gov,  last  accessed  on  May  25,  2005). 


In  the  northern  part  (off  the  Shelter  Island  -  SI/1,  number  1  in  Figure  7),  the  range 
is  between  14°C  and  19.5°C  during  the  year  and  17.5°C  to  19.5°C  during  the  summer.  At 
the  Laurel  Street  Anchorage  (LSA/2,  number  2  in  Figure  7),  it  reaches  21°C,  with  a  range 
of  14°C  to  21°C  during  the  year  but  only  19°C  to  21°C  during  the  summer.  Further  south 
at  the  Bay  Bridge  Anchorage  (BBA/3,  number  3  in  Figure  7)  the  range  is  from  15°C  to 
22°C  and  again  in  the  summer  is  from  20°C  to  22°C.  At  the  Sweet  Water  Channel 
(SWC/4,  number  4  in  Figure  7)  the  range  is  from  14°C  to  25°C  with  22°C  to  25°C  in  the 
summer.  In  Chula  Vista  Marina  (CVM/5,  number  5  in  Figure  7),  it  reaches  a  maximum  of 
26°C  and  a  minimum  of  23°C  in  the  summer.  As  expected,  the  north  part  of  the  bay  is 
colder  than  the  south  and  in  the  summer  the  water  is  much  warmer  than  in  the  winter.  For 
this  study,  a  standard  water  temperature  of  23°C  is  used  since  the  data  were  gathered  in 
the  summer. 
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The  salinity  varies  from  32.5  to  37.5  ppt  and  the  average  is  35  ppt.  The  zonal 
separation  is  similar  to  the  temperature.  In  LSA  and  SI,  the  range  is  from  33  to  35.5  ppt. 
In  BBA,  it  reaches  36  ppt.  In  SWC,  the  range  is  from  32  to  37  ppt.  In  CVM,  it  is  from  34 
to  37.5  ppt.  Again,  detailed  analysis  shows  that  in  the  summer,  salinity  increases  with  an 
average  of  1  ppt  because  of  the  zero  precipitation.  Therefore,  36  ppt  is  used  as  the 
standard  salinity  during  the  summertime. 

E.  ATMOSPHERIC  FORCING 

The  winds  have  a  very  small  effect  on  the  currents  because  of  their  intensity  and 
the  geography  of  the  bay.  Both  mean  westerly  winds  in  the  afternoon  and  mean  easterly 
winds  in  the  morning  and  evening  are  less  than  5  m/s  (10  knots).3  NOAA’s  weather 
description  for  San  Diego  during  June,  July  and  August  points  out  that,  during  these 
months,  there  are  practically  no  storms.  Wind  forcing  is  always  less  significant  than  tidal 
forcing.  The  shallow  waters  make  instrument  deployment  problematic.  With  small 
currents,  the  angular  momentum  of  the  instruments  induces  direction  errors.  Hence,  the 
model  is  more  difficult  to  validate  in  these  areas. 

Annual  precipitation  is  about  0.26  m  (Woodward  and  Clyde,  1996)  and  occurs 
mostly  in  winter.  Therefore,  in  terms  of  estuarine  classification,  San  Diego  Bay  is 
generally  positive,  i.e.,  drainage  inflow  exceeds  evaporation  (Pritchard,  1952).  However, 
during  the  summer,  the  evaporation  rate,  about  0.16  m,  exceeds  precipitation4  (Peeling, 
1975)  and  a  “reversed  estuary”  phenomenon  is  observed  (Defant,  1961).  In  general,  the 
low  inflow  of  fresh  water  provokes  very  small  buoyancy  forcing;  hence  the  density- 
driven  circulation  is  driven  by  seasonal  heating  and  evaporation.  Note  that  for 
precipitation,  in  June,  the  rainfall  is  negligible  averaging  only  0.17  cm,  in  July  0.5  cm  and 
in  August  0.25  cm.5  Therefore,  small  surface  water  mass  flux  (mostly  in  winter)  and 
wind  forcing  for  the  San  Diego  Bay  are  ignored.  This  study  area  is  a  small  basin  with 
circulation  driven  by  tidal  flow  (Fagherazzi  et  ah,  2003). 

3  In  accordance  with  the  Golden  Date  Weather  Services  archive,  the  average  speed  is  between  6.4  and 
7.8  knots.  See  http://ggweather.com,  last  accessed  on  May  25,  2005. 

4  The  precipitation  during  the  summer  is  usually  zero. 

5  NOAA.  See  http://www.wrh.noaa.gov/sgx/climate/san-san-month.htm,  last  accessed  on  May  25, 
2005. 
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The  Californian  Current  has  a  significant  influence  on  the  regional  oceanography 
of  the  Californian  coast.  Under  certain  conditions,  the  alongshore  flow  can  influence  the 
shelf  circulation  in  southern  California,  which  in  turn,  drives  flow  around  Point  Loma. 
This  can  lead  to  either  a  re-circulating  eddy  in  the  lee  of  the  headland  or  northward 
currents  inshore. 

F.  WATER  QUALITY  MONITORING 

In  1960,  an  earthquake  with  a  Richter  scale  of  9  in  Chile  caused  the  biggest 
sudden  rise  in  sea  level  ever  recorded  in  the  San  Diego  area  of  1 .07  m  at  the  Scripps  pier. 
There  is  a  natural  protection  due  to  the  160  km  wide  continental  shelf  of  San  Diego. 
There  is  a  fault  off  San  Diego  Bay,  but  it  is  inactive.  These  are  the  reasons  why  from  the 
15  locally  generated  tsunamis  in  California  since  1812,  only  two  have  occurred  in 
Southern  California,  and  only  one  in  San  Diego,  dating  back  to  1862. 
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Figure  7.  Bay  Wide  Water  Quality  Monitoring  Program  Stations  (From: 

http://www.portofsandiego.org/sandiego  environment/bay  water  sampling. asp), 

last  accessed  on  May  25,  2005. 

Several  mineral  resource  extraction  activities  are  now  occurring  or  have  occurred 
in  the  recent  past  in  San  Diego  Bay.  Among  them,  the  most  important  are  the  production 
of  sodium  chloride  (salt),  bromine  and  other  chemicals  from  sea  water,  magnesium, 
magnesium  compounds  and  brine.  There  is  widespread  toxicity  in  San  Diego  Bay 
sediments  attributable  to  copper,  zinc,  mercury,  polycyclic  aromatic  hydrocarbons, 
polychlorinated  biphenyls  (PCBs)  and  chlordane. 

No  single  chemical  or  chemical  group  has  a  dominant  role  in  contributing  to  the 
identified  toxicity.  Contributions  of  trace  metals  from  vessel  activities  have  long  been 
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suspected  as  a  potentially  large  source  to  San  Diego  Bay.  Shelter  Island  Yacht  Basin,  a 
semi-enclosed  boat  harbor,  has  been  added  to  the  State's  list  of  impaired  water  bodies  (the 
303d  list).  These  contributions  arise  from  specially  formulated  paints,  impregnated  with 
biocides,  and  applied  to  boat  hulls  to  retard  the  growth  of  fouling  organisms  such  as 
barnacles.  Levels  of  potentially  toxic  trace  metals,  especially  copper  and  zinc  (Figure  8) 
exceed  regulatory  water  quality  objectives  and  gradients  of  toxicity  to  mussels  increase  as 
one  enters  the  Shelter  Island  (Figure  1). 
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Figure  8.  Zinc  concentration  in  San  Diego  Bay  on  July  4,  1993  (From: 
http://sdbay.sdsc.edu/html,  last  accessed  on  May  25,  2005). 
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III.  BAROTROPIC  AND  BAROCLINIC  CURRENTS 


Three  Acoustic  Doppler  Current  Profilers  (ADCPs)  (Figure  9)  were  deployed  by 
SPAWAR  in  1993  and  their  measurements  will  be  used  in  this  study.  The  first  was  a 
broad  band  one  (bb)  and  was  deployed  in  position  (32°42’25.8”N,  1 17°13’30.6”  W)  from 
June  22  until  July  23,  1993.  The  second  was  a  narrowband  one  (nbl)  and  was  deployed  in 
position  (32°42’43.98”N,  1 17°12’55.68”W)  from  June  22  until  August  26,  1993.  The  last 
was  a  narrowband  one  (nb2)  and  was  deployed  in  position  (32°  42’  17.22”N,  117°10’ 
8.88”W)  from  June  23  until  August  27,  1993. 


Figure  9.  Positions  of  the  ADCPs  deployed  in  1993  (From:  WQMAP). 

By  checking  the  ADCP  data  inside  the  bay  (nbl,  nb2),  San  Diego  Bay  was 
concluded  to  be  vertically  well  mixed.  The  u  and  v  components  are  plotted  in  three 
different  depths  for  the  ADCP  measurements  to  prove  the  aforementioned  statement 
(Figures  10-15). 
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V  (cm/s)  3  U  (cm/s) 


U  in  water  column  for  NB1 


time  (in  10  s) 

U  component  from  ADCP  data  (nbl)  for  surface  (yellow),  middle  depth 
(purple)  and  bottom  (blue). 


V  in  water  column  for  NB1 


time  (in  10  s) 

Figure  11.  V  component  from  ADCP  data  (nbl)  for  surface  (yellow),  middle  depth 

(purple)  and  bottom  (blue). 
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U  in  water  column  for  NB2 


time  (in  10  s) 

Figure  12.  U  component  from  ADCP  data  (nb2)  for  surface  (yellow),  middle  depth 

(purple)  and  bottom  (blue). 


V  in  water  column  for  NB2 


time  (in  10  s) 


Figure  13.  V  component  from  ADCP  data  (nb2)  for  surface  (yellow),  middle  depth 

(purple)  and  bottom  (blue). 
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u  (cm/s)  3  V  (cm/s) 


V  in  water  column  for  BB 


time  (in  10  s) 

V  component  from  ADCP  data  (bb)  for  surface  (yellow),  middle  depth 
(purple)  and  bottom  (blue). 


U  in  water  column  for  BB 


time  (in  10  s) 

Figure  15.  U  component  from  ADCP  data  (nb2)  for  surface  (yellow),  middle  depth 

(purple)  and  bottom  (blue). 
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For  nbl,  the  correlation  coefficient  between  the  surface  and  bottom  is  97.16%  for 
the  u-component  and  97.16%  for  the  v-component  96.32%.  For  nb2,  the  correlation 
coefficient  between  the  surface  and  bottom  is  91.89%  for  the  u-component  and  94.71% 
for  the  v-component  96.32%.  Between  the  surface  and  bottom,  the  phase  matches  very 
well  and  so  does  the  trend  in  the  respective  amplitudes.  In  terms  of  actual  amplitude 
matching,  it  is  possible  to  optically  verify  the  good  match.  The  differences  observed 
cannot  contradict  that  San  Diego  Bay  is  well  mixed. 

The  bb  is  a  more  interesting  case  because  it  is  not  in  the  main  bay,  since  it  is 
further  south  from  Coronado  Island  to  the  east  of  Point  Loma.  In  that  case,  an  excellent 
correlation  still  exists  for  the  v  component  but  not  for  the  u.  The  correlation  coefficient 
between  the  surface  and  bottom  is  92.94%  for  the  v-component  and  35.19%  for  the  u- 
component.  By  checking  the  correlation  between  the  u  component  in  the  surface  and  in 
middle  depth,  the  coefficient  improves  and  becomes  78.44%  and  by  comparing  the  v 
component  of  the  bottom  and  the  middle  depth  it  becomes  61.65%.  In  order  to  understand 
the  reasons  for  this,  the  power  spectrum  diagram  (PSD)  of  u  component  at  bb  is 
calculated  for  both  the  surface  and  the  bottom  (Figures  16-17). 


PSD  Diagram  fof  bb  U  component  surface 


Figure  16.  Power  spectrum  diagram  (PSD)  for  bb  u  component  in  the  surface. 
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PSD  Diagram  for  bb  U  component  bottom 


Figure  17.  Power  spectrum  diagram  (PSD)  for  bb  u  component  at  the  bottom. 

Distinction  of  the  power  spectra  of  u  at  bb  between  the  surface  and  bottom  clearly 
indicate  that  waters  cannot  be  treated  as  well-mixed.  Obviously,  the  spikes  both  at  low 
and  high  frequencies  (lower  than  1  cycle  per  day  and  higher  than  2  -  especially  at  3  and  4 
cycles  per  day)  are  very  significant.  Since  it  is  not  in  a  more  “open  ocean”  environment, 
it  is  very  normal  that  both  the  California  Current  System  and  the  wind  play  a  significant 
role.  However,  it  is  not  only  the  long  tenn  or  short  term  forcing  that  is  responsible  for  this 
difference.  Even  when  filtering  the  data  to  keep  only  the  part  that  results  from  the  tidal 
influence,  the  correlation  was  not  satisfactory.  A  band-pass  butterworth  filter  was  used  to 
filter  out  all  frequencies  either  smaller  than  0.5  cycles  per  day  or  greater  than  2.5  cycles 
per  day.  For  v,  an  absolutely  perfect  correlation  (100%)  was  achieved,  but  for  u,  a 
correlation  coefficient  of  only  49.53%  was  achieved. 

Since  the  correlation  of  the  v  component  is  perfect,  something  obviously 
influences  the  u  component.  It  is  possible  to  easily  conclude  the  proximity  of  the  ADCP 
data  to  the  Zuniga  jetty  (which  has  a  North-South  direction)  influences  the  u  component 
due  to  reflections.  Hence,  suffice  it  to  say  that  there  are  well  mixed  waters  inside  San 
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Diego  Bay  and  without  ruling  out  the  well-mixed  nature  of  the  waters  in  the  mouth  of  the 
Bay,  it  is  possible  to  conclude  that  the  waters  are  less  well-mixed.  Hence,  the  application 
of  a  3-D,  if  not  necessary,  will  produce  better  results.  However,  firstly,  it  was  concluded 
that  the  influence  of  other  factors  (such  as  the  wind  driven  circulation)  cannot  be  entirely 
ignored  outside  the  bay,  and  secondly,  that  there  are  well  mixed  waters  only  inside  San 
Diego  Bay. 


ADCP  site 

Velocity 

component 

Correlation 

coefficient 

Remarks 

Nbl 

u 

97.16% 

unfiltered 

Nbl 

V 

96.32% 

unfiltered 

Nb2 

u 

91.89% 

unfiltered 

Nb2 

V 

94.71% 

unfiltered 

bb 

u 

35.19% 

unfiltered 

bb 

V 

92.94% 

unfiltered 

bb 

u 

49.53% 

Filtered  0.5  -  2.5  cycles  per 
day 

bb 

V 

100.00% 

Filtered  0.5  -  2.5  cycles  per 
day 

Table  2.  Correlation  between  surface  and  bottom  measurements  in  all  ADCP  sites  - 

proof  of  barotropic  nature  of  San  Diego  Bay. 
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IV.  HYDRODYNAMIC  MODEL 


A.  MODEL  DESCRIPTION 

The  hydrodynamic  model  used  in  this  study  is  the  “Water  Quality  Management 
and  Analysis  Package”  (WQMAP)  developed  by  the  Applied  Science  Associates,  Inc.  in 
Narragansett,  Rhode  Island.  It  consists  of  a  boundary  confonning  grid  generation  model, 
a  three-dimensional6  time  dependent  primitive  equation  model,  and  a  suite  of  pollutant 
transport  and  fate  models.  These  models  are  executed  on  a  boundary  fitted  grid  system. 
However,  they  can  be  operated  on  any  orthogonal  curvilinear  grid  or  a  rectangular  grid, 
which  are  special  cases  of  the  boundary  fitted  grid.  The  model  is  configured  to  run  in  a 
vertically  averaged  (barotropic)  mode  or  in  a  fully  three-dimensional  mode.  Several 
approximations  are  made  in  the  model  fonnulation,  including  the  hydrostatic  and 
Boussinesq  approximations  (Chu  et  ah,  2004).  As  happens  in  all  vertically-averaged 
models,  it  is  implicitly  assumed  that  the  water  density  and  velocity  are  nearly  depth- 
independent  over  the  water  column.7  It  is  also  assumed  that  the  water  depth  is  sufficiently 
shallow  so  as  to  neglect  the  Ekman  layer  (Cheng  et  al.,  1983). 

WQMAP  embeds  the  Geographic  Information  System  (GIS)  with  maps  for  any 
region  in  the  world  and  environmental  data  management  tools.  The  GIS  can  be  used  to 
provide  input  data  to  the  models,  as  well  as  to  present  model  predictions  better.  The 
model  solves  the  three  dimensional  conservation  of  mass,  momentum,  and  energy 
equations  on  a  spherical,  non-orthogonal,  boundary  confonning  grid  system  and  can  be 
applied  to  both  estuarine  and  littoral  regions.  The  eddy  viscosities  are  either  specified  by 
the  user  or  based  on  a  one-equation  turbulent  kinetic  energy  (TKE)  model.  Output  of  the 
TKE  model  can  be  used  in  conjunction  with  a  prescribed  mixing  length  to  determine  the 
vertical  eddy  viscosities.  The  model  predicts  time  varying  fields  of  surface  elevations  and 
velocities.  Environmental  forcing  inputs  include  tides,  river  inflows,  surface  elevations, 

6  WQMAP  can  also  run  in  a  two-dimensional  mode,  which  is  a  better  simplifying  choice  for  San 
Diego  Bay. 

7  Normally,  the  horizontal  density  gradients  are  treated  explicitly  in  the  momentum  equations; 
however,  due  to  the  predominantly  tidal  current  production  in  San  Diego  Bay,  the  horizontal  density 
gradient  can  be  neglected  in  short-time  prediction. 
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and  wind  fields.  A  sigma  stretching  system  is  incorporated  to  diagram  the  free  surface 
and  bottom  to  resolve  bathymetric  variations  (Annstrong,  2004).  Calculations  are 
achieved  on  a  space-staggered  grid  system  in  the  horizontal  and  a  non-staggered  system 
in  the  vertical. 

The  three-dimensional  hydrodynamic  equations  contain  fast  moving  external 
gravity  waves  and  slow  moving  internal  gravity  waves.  The  equations  of  motion  are  split 
into  vertically  averaged  equations  (exterior  mode)  and  vertical  structure  equations 
(interior  mode).  This  technique  allows  the  calculation  of  the  free  surface  elevation  from 
the  exterior  mode  and  the  three  dimensional  currents  and  thermodynamic  properties  from 
the  interior  mode.  The  external  mode  equations  are  obtained  by  integrating  the  three 
dimensional  equations  of  continuity  and  conservation  of  momentum  from  a  =  0  to  a  =  1 
(WQMAP  Technical  Manual,  2003).  In  the  exterior  mode,  the  Helmholtz  equation,  given 
in  terms  of  the  sea  surface  elevation,  is  solved  to  ease  the  time  step  restrictions  normally 
imposed  by  the  aforementioned  gravity  wave  propagation.  In  the  interior  mode,  the  flow 
is  predicted  by  an  explicit  finite  difference  method,  except  that  the  vertical  diffusion  term 
is  treated  implicitly  (Madala  and  Piaczek,  1977).  The  time  step  generally  remains  the 
same  for  both  exterior  and  interior  modes  (Spaulding  et  ah,  1999).  Hence,  the  user  can 
specify  the  variable  (currents,  temperature,  surface  elevation  and  salinity)  over  a  certain 
simulation  period  and  then  the  results/  predictions  are  presented  in  either  color  contours 
or  vectors. 

B.  MOMENTUM  AND  CONTINUITY  EQUATIONS  FOR  BAROTROPIC 

MODE 

If  U  is  the  vertically  averaged  component  in  £,  -direction,  V  the  vertically 
averaged  component  in  r\  -direction,  Ah  the  horizontal  eddy  viscosity,  xw  the  wind 

stress,  and  xh  the  bottom  stress,  then  the  two-dimensional  vertically  averaged 
conservation  of  momentum  equations  are  as  follows: 
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The  continuity  equation  is  given  by 
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Here,  pG  is  the  characteristic  water  density  (1,025  kg/m'),  /  is  the  Coriolis 
parameter  (f=  2Qsin(p),  R  is  the  earth  radius,  and 

z  +  H 
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t  +  H 

where  H  is  the  water  depth,  D  =  H  +  C  ,  is  the  total  water  depth  (bathymetry  and  surface 
elevation),  a  =  1  represents  the  ocean  surface,  and  a  =  0  represents  the  ocean  bottom. 
Finally,  the  transport  equation  for  \//  (temperature,  salt  or  any  substance)  is  as  follows, 
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C.  COORDINATE  SYSTEM  AND  BOUNDARY  CONDITIONS 

An  important  feature  of  WQMAP  is  its  hybrid  orthogonal  curvilinear-terrain 

following  coordinate  system.  Let  (cp,0)  be  the  latitude  and  longitude  and  (c,,il)  be  the 

coordinates  on  a  generalized  orthogonal  curvilinear  coordinate  system.  Let  C,  be  the 
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surface  elevation  and  a  the  vertical  o-coordinate.  The  two  metric  tensors  connecting 
((p,0)to  (c,r|)  are  defined  and  the  coefficient  gn  is  the  metric  tensor  in  c  -direction  and 


the  coefficient  g22  metric  tensor  in  q  -direction.  These  tensors  permit  the  model  to 

transfonn  the  boundary  fitted  grid  to  a  numerical  grid  employed  for  spatial  discretization 
utilized  in  an  Arakawa  C  Grid  (Chu  et  al.,  2004) 


gn 


(6) 


which  are  the  metric  tensors  in  £,  and  r)  directions.  Furthermore,  at  the  surface, 
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and  at  the  bottom, 
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The  land  boundaries  are  assumed  impermeable.  Therefore,  the  normal  component 
of  the  velocity  is  set  to  zero.  Sea  surface  elevation  or  tidal  harmonic  constituents  can  be 
specified  as  a  function  of  time  along  the  open  boundaries.  At  the  closed  boundaries,  the 
transport  of  substance  (i.e.,  salinity)  is  zero.  At  the  open  boundaries,  the  concentration  is 
specified  during  the  inflow,  using  the  characteristic  values  for  San  Diego  Bay  (for 
temperature,  23°  for  summer  months  is  used  and  for  salinity  36  ppt  is  used).  On  outflow, 
the  substance  is  advected  out  of  the  model  domain  according  to  the  following  equations: 
8¥  .  u  dy/  dy/  v  d y/ 
dt  ’  8t  dl) 

D.  RECTANGULAR  GRID 

The  scenario  used  for  this  research  uses  both  the  two-dimensional  and  three- 
dimensional  hydrodynamic  model,  taking  into  account  the  minimal  wind  considerations 
for  San  Diego  Bay,  as  analyzed  in  Chapter  I.  WQMAP  is  supplied  with  maps  of  the  area. 
A  small  spatial  extent  ASA  file  (.bdm),  containing  high  resolution  characteristics,  is  used. 
The  grid  used  (Figure  18)  was  created  with  a  detailed  resolution  (Imax  and  Jmax)  150  by 
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200  (30,000  nodes)  grid.  For  this  grid,  the  average  horizontal  resolution  is  40m.  The 
rectangular  grid  used  improves  the  previously  created  boundary-fitted  ones  (Armstrong, 
2004)  having  a  better  resolution  and  realistically  takes  into  consideration  the  Zuniga  jetty 
(Figure  19),  not  included  in  previous  studies. 
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Figure  19.  Grid  representation  taking  into  account  the  Zuniga  jetty  (From:  WQMAP). 

The  bathymetry  is  from  the  NOAA  3-second  data  base  supported  by  the  National 
Ocean  Service  (NOS)  soundings  database  (i.e.,  Hydrographic  Survey  Data,  version  3.2) 
averaged  within  each  cell.  Hence,  it  is  adequate  for  model  resolution;  however,  it  should 
be  applied  with  caution  near  the  coastline  because  the  resolution  of  the  data  base  does  not 
conform  to  small  details  in  the  coastline.  As  seen  in  Chapter  V,  this  is  a  factor  to  be  taken 
into  account  since  the  current  meters  are  located  very  close  to  the  shore.  It  should  also  be 
mentioned  that  the  model  requests  a  minimum  depth  (set  at  2  m  in  our  study).  There  is 
also  a  smoothing  of  the  bathymetry,  which  is  essential  for  the  model  function,  but  might 
cause  significant  errors  in  very  shallow  areas. 
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V.  MODEL  EVALUATION 


The  tidal  forcing  in  the  entrance  of  the  bay  is  reproduced  with  time  series  from 
Wxtide32  model  database  based  on  verified  data  from  the  NOOA/NOS  CO-OPS.  All 
measurements  (hourly  height  water  level)  had  as  a  reference  the  MLLW  and  were 
conducted  between  June  22,  1993  and  August  27,  1993.  Figure  20  shows  the  water 
elevation,  in  accordance  with  NOAA  San  Diego  Station  number  9410170,  located  at 
(32°42’48”N,  1 17°10’24”W).. 

The  aforementioned  tidal  forcing  (Table  3)  was  chosen  to  verify  the  model  with 
actual  data  available  for  the  certain  period  of  time. 

The  positions  of  the  ADCPs  deployed  will  be  repeated:  (1)  the  broad  band  ADCP 
(at  bb)  at  (32°42’25.8”N,  117°13’30.6”  W)  from  June  22  until  July  23,  1993,  (2) 
narrowband  ADCP  (at  nbl)  at  (32°42’43.98”N,  1 17°12’55.68”W)  from  June  22  until 
August  26,  1993,  and  (3)  the  narrowband  ADCP  (at  nb2)  at  (32°42’  17.22”N,  117°10’ 
8.88”W)  from  June  23  until  August  27,  1993.  An  initial  comparison  of  unfiltered  data 
and  original  model  run  is  conducted  to  see  if  they  are  correlated.  Therefore,  the  average 
U  and  V  components  of  the  ACDP  data  and  the  initial  model  results  are  plotted  (Figures 
21-24). 
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MM 

0.0015122 

0.0056 

0.009 

243.30 

129.84 

0.4 

MSF 

0.0028219 

0.0036 

0.010 

249.32 

178.00 

0.14 

ALP1 

0.0343966 

0.013 

0.0024 

330.95 

231.81 

0.034 

2Q1 

0.0357064 

0.0027 

0.012 

202.55 

233.41 

0.055 

*Q1 

0.0372185 

0.0345 

0.016 

67.28 

29.68 

4.8 

*01 

0.0387307 
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0.017 

307.64 
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1.3e+002 

NOl 
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0.012 
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0.23 
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001 
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0.0047 

0.013 
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EPS2 
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224.64 

229.71 
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*MU2 
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3.9e+002 

ETA2 
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0.0008 

0.007 
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239.26 

0.013 

*M03 

0.1192421 
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M3 

0.1207671 

0.0005 
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243.70 
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1.9 

*MK3 

0.1222921 

0.0019 
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90.76 

8.88 
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SK3 

0.1251141 

0.0004 

0.000 

158.31 

49.29 

1.8 

*MN4 

0.1595106 

0.0014 

0.001 

230.16 

20.69 

7.6 

*M4 

0.1610228 

0.0044 

0.001 

113.75 

7.48 

67 

SN4 

0.1623326 

0.0006 

0.001 

214.31 

50.35 

1.5 

*MS4 

0.1638447 

0.0032 

0.000 

144.46 

10.24 

41 

S4 

0.1666667 

0.0004 

0.000 

200.58 

87.45 

0.85 

*2MK5 

0.2028035 

0.0006 

0.000 

311.82 

30.51 

5.3 

2SK5 

0.2084474 

0.0002 

0.000 

270.56 

76.50 

0.68 

*2MN6 

0.2400221 

0.0002 

0.000 

80.78 

25.22 

6 

M6 

0.2415342 

0.0001 

0.000 

9.97 

75.96 

0.84 

*2MS6 

0.2443561 

0.0001 

0.000 

300.71 

37.62 

2.9 

*2SM6 

0.2471781 

0.0001 

0.000 

288.67 

33.45 

2.8 

*3MK7 

0.2833149 

0.0002 

0.000 

319.58 

26.45 

5.2 

M8 
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0.0001 

0.000 

5.05 

57.00 

1.4 

Table  3.  Harmonic  decomposition  of  initial  forcing. 
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Figure  20.  Tidal  forcing  used  for  the  verification  of  the  model  (From:  http://www.co- 

ops.nos.noaa.gov,  last  accessed  on  May  25,  2005). 
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Figure  2 1 .  Average  U  component  from  ADCP  (nbl)  data  (purple  line)  compared  to 

initial  model  results  (blue  line). 
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V  for  nbl 


Time  (date) 


Figure  22.  Average  V  component  from  ADCP  (nbl)  data  (purple  line)  compared  to 

initial  model  results  (blue  line). 

U  for  nb2 


Time  (date) 


Figure  23.  Average  U  component  from  ADCP  (nb2)  data  (purple  line)  compared  to 

initial  model  results  (blue  line). 
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V  for  nb2 


Time  (date) 

Figure  24.  Average  V  component  from  ADCP  (nb2)  data  (purple  line)  compared  to 

initial  model  results  (blue  line). 

Statistical  analysis  shows  that  a  good  correlation  exists  between  most  of  the 
compared  signals  in  both  phase  and  amplitude  (correlation  coefficient  above  90%  in  all 
cases).  For  nbl,  the  u  speed  between  the  data  and  the  model  has  a  correlation  coefficient 
of  91.87%  and  can  be  verified  optically.  The  observational  u-velocity  ranges  between  - 
51.8  and  44.5  cm/s  and  the  modeled  u-velocity  changes  between  -46.9  and  40.8  cm/s. 
The  difference  between  the  observational  and  modeled  mean  u-velocity  is  0.49  cm/s. 
Furthermore,  the  root  mean  square  error  (between  model  and  observation)  is  9.02  cm/s. 
For  nbl,  the  v  speed  between  the  data  and  the  model  has  a  correlation  coefficient  of 
91.66%  and  can  be  verified  optically.  The  observational  v-velocity  ranges  between  -31.6 
and  29.6  cm/s  and  the  modeled  v-velocity  changes  between  -37.0  and  32.0  cm/s.  The 
difference  between  the  observational  and  modeled  mean  v-velocity  is  -0.65  cm/s.  The 
root  mean  square  error  of  v-velocity  is  6.83  cm/s. 

For  nb2,  the  u  speed  between  data  and  model  has  a  correlation  coefficient  of 
92.60%  and  can  be  verified  optically.  The  observational  u-velocity  ranges  between  -42.8 
and  32.8  cm/s  and  the  modeled  u-velocity  changes  between  -40.2  and  34.3  cm/s.  The 
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difference  between  the  observational  and  modeled  mean  u-velocity  is  1.0862  cm/s. 
Furthermore,  the  root  mean  square  error  (between  model  and  observation)  is  6.7356  cm/s. 
For  nbl,  the  v  speed  between  the  data  and  the  model  has  a  correlation  coefficient  of  92.60 
%  and  can  be  verified  optically.  The  observational  v-velocity  ranges  between  -34  and 
42.7  cm/s  and  the  modeled  v-velocity  changes  between  -20.4  and  23.5  cm/s.  The 
difference  between  the  observational  and  modeled  mean  v-velocity  is  -1.2971.  The  root 
mean  square  error  of  v-velocity  is  8.5035.  In  order  to  verify  the  minimal  significance  of 
the  wind  driven  circulation,  the  power  spectra  diagrams  (PSD)  of  u  and  v  for  nbl  and  nb2 
are  presented  (Figures  25  -  28). 

x  105  PSD  diagram  for  nbl  U  compnent 
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Figure  25.  Power  Spectrum  Diagram  for  nbl  u  component. 
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PSD  Diagram  for  nbl  V  component 


x  105 


PSD  Diagram  for  nb2  U  component 


Figure  27.  Power  Spectrum  Diagram  for  nb2  u  component. 
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PSD  Diagram  for  nb2  V  component 
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Figure  28.  Power  Spectrum  Diagram  for  nb2  v  component. 

It  is  obvious  that  the  long-term  wind  driven  effect  is  minimal  and  can  be 
neglected.  The  circulation  is  due  mainly  to  the  semi-diurnal  and  diurnal  tides,  as 
expected.  Figures  29  through  32  are  identical  during  the  integration  period  for  both  nbl 
and  nb2  with  peaks  at  1  cycle  (diurnal  tides)  and  2  cycles  per  day  (semi-diurnal  tides). 
Long  term  effects  (less  than  1  cycle  per  day)  or  short  tenn  ones  (more  than  half  a  cycle 
per  day)  cannot  be  identified  since  there  is  very  low  energy  in  any  other  frequency. 

However,  it  will  be  statistically  proven  that  wind  is  insignificant.  For  this  reason, 
three  separate  model  runs  will  be  used.  The  first  will  be  with  tidal  forcing  and  an  average 
westerly  wind  of  5m/s;  the  second  with  only  tidal  forcing;  and  the  third  only  with  wind. 
The  comparisons  between  the  runs,  which  include  the  tidal  forcing,  are  practically 
identical  irrespective  of  the  wind  influence.  For  this  reason,  it  is  not  necessary  to  present 
all  the  results.  Note  that  for  nbl  u  component  the  same  maximum  and  minimum  value 
exists  in  both  cases.  The  correlation  coefficient  is  100%,  the  difference  of  the  mean 
values  is  -0.0037  cm/s  and  the  difference  between  the  standard  deviations  is  -0.0025.  The 
root  mean  square  error  is  as  low  as  0.0649  cm/s. 


0.5  1  1.5  2  ZtT  3  3  5  4  4iT  5 

Frequency  (Cycles  per  day) 
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The  comparison  as  regards  the  run  that  takes  into  account  the  wind  without  tidal 
forcing  with  the  one  that  includes  tidal  forcing  is  very  clear.  Again  if  we  take  as  an 
example  the  nbl  u  component,  the  maximum  value  with  all  forcing  is  40.8  cm/s  and  only 
with  wind  is  0.7  cm/s  and  the  minimum  -46.9  cm/s  and  -0.1  m/s  respectively.  There  is  no 
correlation  at  all  (correlation  coefficient  0).  The  contribution  of  the  wind  to  the  current 
speed  can  be  plotted.  Actually,  in  Figure  29,  the  u  component  is  plotted  for  ADCP  nbl  in 
the  run  with  wind  only  and  in  the  run  with  total  forcing.  Flence,  we  can  verify  again  that 
wind  is  insignificant  in  San  Diego  Bay  forcing  and  therefore  by  neglecting  it,  our  results 
are  entirely  valid.  Tables  4-9  provide  detailed  statistical  analysis  on  model  performance. 


Wind  contribution  to  total  forcing 


Time  (date) 

Figure  29.  Average  u  component  from  ADCP  (nbl)  only  with  wind  forcing  (purple 
line)  compared  to  model  results  in  a  total  forcing  run  (blue  line). 
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DATA 

MODEL 

DIFFERENCE 

COMPARISON 

max 

44.5 

40.8 

-3.7 

min 

-51.8 

-46.9 

4.9 

mean 

-0.5384 

-0.0483 

0.4901 

std 

22.7892 

21.3238 

-1.4654 

nne 

-0.9103 

nnse 

9.0238 

ecv 

-16.7615 

corrl 

91.87% 

Table  4.  Comparison  of  u  component  for  nb  1 . 


DATA 

MODEL 

DIFFERENCE 

COMPARISON 

max 

29.60 

32.00 

2.40 

min 

-31.60 

-37.00 

-5.40 

mean 

0.6190 

-0.0315 

-0.6505 

std 

14.2799 

16.7878 

2.5079 

nne 

1.0509 

nnse 

6.8315 

ecv 

11.0361 

coni 

91.66% 

Table  5 .  Comparison  of  v  component  for  nb  1 . 


DATA 

MODEL 

DIFFERENCE 

COMPARISON 

max 

60.51 

59.74 

-0.77 

rnin 

0.10 

0.00 

-0.10 

mean 

23.1499 

23.4548 

0.3049 

std 

13.7075 

13.6486 

-0.0589 

nne 

0.0132 

rmse 

9.7512 

ecv 

0.4212 

coni 

74.61% 

Table  6.  Comparison  of  speed  for  nb  1 . 


DATA 

MODEL 

DIFFERENCE 

COMPARISON 

max 

32.80 

34.30 

1.50 

rnin 

-42.80 

-40.20 

2.60 

mean 

-0.9540 

0.1322 

1.0862 

std 

17.0167 

17.4624 

0.4457 

nne 

-1.1386 

nnse 

6.7356 

ecv 

-7.0605 

coni 

92.60% 

Table  7.  Comparison  of  u  component  for  nb2. 
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DATA 

MODEL 

DIFFERENCE 

COMPARISON 

max 

42.70 

23.50 

-19.20 

min 

-34.00 

-20.40 

13.60 

mean 

1.2354 

-0.0617 

-1.2971 

std 

16.9599 

10.2599 

-6.70 

line 

1.05 

nnse 

8.5035 

ecv 

6.8833 

corrl 

92.60% 

Table  8.  Comparison  of  v  component  for  nb2. 


DATA 

MODEL 

DIFFERENCE 

COMPARISON 

max 

60.32 

46.56 

-13.76 

rnin 

0.10 

0.00 

-0.10 

mean 

20.2309 

17.4891 

-2.7418 

std 

12.7761 

10.5454 

-2.2307 

line 

0.1355 

rmse 

8.0369 

ecv 

0.3973 

corrl 

80.66% 

Table  9.  Comparison  of  speed  for  nb2. 


Comparison  of  the  current  speed  between  observation  and  the  model  is 
acceptable.  As  regards  the  point  south  of  Coronado  Island  (bb),  when  fdtering  the  data  to 
keep  only  the  part  that  results  from  the  tidal  influence,  an  improved  correlation  between 
the  model  and  observation  occur.  Again,  a  band-pass  butterworth  filter  is  first  applied  to 
filter  out  all  frequencies  either  smaller  than  0.5  cycles  per  day  or  greater  than  2.5  cycles 
per  day.  The  correlation  coefficient  is  85.09%  for  u-velocity  and  78.49%  for  v-velocity. 
Optically,  the  amplitude  difference  between  the  model  and  observation  improves  as  well. 
Thus,  it  is  possible  to  conclude  that  by  filtering  the  results  of  the  non-tidal  influence,  a  2- 
D  model  can  be  used  in  the  entrance  of  San  Diego  Bay  with  relatively  better  results. 

The  middle  and  upper  levels  contribute  much  more  to  the  vertical  average  u 
component  of  bb  than  the  lower  (to  bottom)  depths  and  indicate  its  behavior  much  better. 
Actually  the  correlation  coefficient  of  the  average  u  component  and  the  surface  u  is 
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91.63%  and  between  the  average  u  component  and  the  middle  depth  u  is  even  higher 
(92.14%).  Moreover,  the  middle  depth  represents  the  u  component  better  since  it  has  also 
a  smaller,  relative  mean  error  (RME) 

n 

^  j^rmod  X°bS  ^ 

RME  =  — - 

n 


i= 1 


root  mean  square  error  (RMSE) 


RMSE  = 


i 


^  mod  X°bs  ^  ^ 


and  the  error  coefficient  of  variation  (ecv), 


ECV  = 


RMSE 


where  (X”od,X°bs)  are  the  modeled  and  observational  data. 

The  bottom  u  component  and  the  average  u  have  a  correlation  coefficient  of  only 
62.22%.  For  the  v  component,  the  correlation  between  the  average  and  the  v  component 
in  all  depths  is  high:  For  bottom,  the  correlation  coefficient  is  96.24%,  for  surface 
99.26%,  and  again  for  middle  depth,  it  is  the  best  with  a  coefficient  of  99.43%  and  better 
nne,  rinse  and  ecv  than  the  other  cases.  In  order  to  evaluate  the  model  further,  the  water 
elevation  in  position  nb2  is  calculated  and  despite  the  fact  that  exact  verified  data  is  not 
available,  an  attempt  will  be  made  to  compare  it  in  two  different  ways.  First,  Table  10 
describes  all  results  after  harmonic  decomposition. 

To  compare  the  model  results  with  a  previous  study  (Wang  et  ah,  1998),  the  four 
main  tidal  constituents  (Kl,  01,  S2,  M2)  at  downtown  (further  north  and  deeper  than  the 
nb2  position)  are  used.  Amplitude  difference  (model  minus  observation)  is  3.83  cm  for 
M2  (but  only  1.08  cm,  if  compared  to  that  of  south  San  Diego),  3.73  cm  for  Ki,  -2.19  cm 
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for  Oi,  and  1.1  cm  for  S2.  Phase  difference  is  1.71°  for  M2,  25.94°  for  Ki,  45.33°  for  Oi, 
and  5.41°  for  S2.  Differences  in  bathymetry  and  analysis  can  justify  such  small 
discrepancies,  especially  since  there  is  no  certainty  about  the  position  for  which  the  data 
was  collected. 


TIDE 

FREQ 

AMPL 

AMP. ERR 

PHASE 

PH  ERR 

SNR 

MSF 

0.0028219 

0.0061 

0.014 

266.81 

160.24 

0.2 

*2Q1 

0.0357064 

0.0067 

0.003 

337.90 

21.94 

6.8 

*Q1 

0.0372185 

0.0364 

0.002 

76.81 

4.66 

2.4e+002 

*01 

0.0387307 

0.1952 

0.003 

125.14 

0.86 

5e+003 

*N01 

0.0402686 

0.0096 

0.003 

19.17 

16.04 

13 

*K1 

0.0417807 

0.3773 

0.002 

60.54 

0.44 

2.5e+004 

J1 

0.0432929 

0.0026 

0.002 

97.99 

69.16 

1.2 

*001 

0.0448308 

0.0157 

0.002 

129.23 

9.23 

40 

*UPS1 

0.0463430 

0.0047 

0.003 

289.53 

32.85 

2.7 

*N2 

0.0789992 

0.1226 

0.014 

203.96 

7.74 

75 

*M2 

0.0805114 

0.5804 

0.015 

270.27 

1.36 

1.6e+003 

*S2 

0.0833333 

0.2144 

0.013 

267.04 

3.93 

2.6e+002 

ETA2 

0.0850736 

0.0077 

0.011 

7.45 

98.73 

0.48 

*M03 

0.1192421 

0.0042 

0.001 

258.54 

22.76 

8.5 

*M3 

0.1207671 

0.0021 

0.001 

172.85 

40.42 

2.4 

*MK3 

0.1222921 

0.0085 

0.001 

219.46 

10.21 

33 

*SK3 

0.1251141 

0.0026 

0.001 

208.56 

32.29 

3.7 

*MN4 

0.1595106 

0.0039 

0.002 

15.82 

21.09 

6.1 

*M4 

0.1610228 

0.001 

75.84 

8.11 

71 

*MS4 

0.1638447 

0.002 

71.22 

11.13 

23 

S4 

0.1666667 

0.0014 

0.001 

66.29 

51.91 

1.2 

*2MK5 

0.2028035 

0.0037 

0.001 

185.30 

24.66 

6.6 

2SK5 

0.2084474 

0.0003 

0.001 

258.37 

225.54 

0.057 

2MN6 

0.2400221 

0.0029 

0.002 

355.66 

52.92 

1.5 

*M6 

0.2415342 

0.0059 

0.002 

52.23 

22.32 

6.5 

*2MS6 

0.2443561 

0.0080 

0.003 

72.33 

18.29 

10 

2SM6 

0.2471781 

0.0019 

0.002 

83.37 

75.50 

0.65 

*3MK7 

0.2833149 

0.0042 

0.002 

108.25 

31.28 

3.4 

*M8 

0.3220456 

0.0007 

0.000 

295.35 

30.49 

3.3 

Table  10.  Harmonic  decomposition  of  modeled  results. 


The  best  choice  for  a  valid  comparison  is  to  compare  the  model  results  at  nb2 
with  the  NOAA  San  Diego  station  at  the  pier.  The  data  is  gathered  too  close  to  land,  and 
therefore,  cannot  be  exactly  compared.  However,  it  can  be  compared  with  a  point  very 
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near  the  actual  San  Diego  Station  position.  Amplitude  difference  (model  minus 
observation)  is  2.51  cm  for  M2,  0.94  cm  for  Ki,  0.84  cm  for  Oi,  and  0.71  cm  for  S2.  Phase 
difference  is  0.75°  for  M2,  26.08°  for  Ki,  29.58°  for  Oi,  and  48.96°  for  S2.  Again,  it  can  be 
stated  that  differences  in  bathymetry  and  analysis  can  justify  such  small  discrepancies, 
especially  since  there  is  no  precision  about  the  position  for  which  the  data  was  compared. 
The  results  from  this  more  accurate  comparison  are  even  better  than  the  previous  one. 

Overall,  the  model  results  are  reasonably  good,  especially  taking  into  account  that 
the  comparison  between  data  and  model  is  not  at  exactly  the  same  position  and  the 
proximity  of  the  ADCPs  to  the  shore.  If  finer  grid  and  more  accurate  bathymetry  are 
used,  the  model  results  may  be  further  improved. 


NOAA  DATA 
(more  precise) 

SPAWAR  DATA 

M2  (amplitude  difference) 

+  2.51  cm 

+  3.83  cm 

Ki  (amplitude  difference) 

-  0.94  cm 

+  3.73  cm 

0 1  (amplitude  difference) 

-  0.84  cm 

-  2.19  cm 

S2  (amplitude  difference) 

+  0.71  cm 

-1.1  cm 

Table  1 1 .  Comparison  of  model  elevation  amplitude  with  NOAA  and  SPAWAR 

results. 


NOAA  DATA 
(more  precise) 

SPAWAR  DATA 

M2  (phase  difference  in  degrees) 

+  0.75 

-1.71 

Ki  (phase  difference  in  degrees) 

-  26.08 

-25.94 

0 1  (phase  difference  in  degrees) 

+29.58 

-45.33 

S2  (phase  difference  in  degrees) 

-48.96 

+  5.41 

Table  12.  Comparison  of  model  elevation  phase  with  NOAA  and  SPAWAR  results. 
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VI.  HYDROCHEMICAL  MODEL 


A.  MODEL  DESCRIPTION 

CHEMMAP  is  a  chemical  spill  model,  which  predicts  the  trajectory  and  fate  of 
floating,  sinking,  evaporating,  soluble  and  insoluble  chemicals  and  product  mixtures.  It 
can  estimate  the  distribution  of  chemical  elements  (as  mass  and  concentrations)  on  the 
surface,  in  the  water  column  and  in  the  sediments.  The  model  is  3D,  separately  tracking 
surface  slicks,  entrained  droplets  or  particles  of  pure  chemical,  chemical  adsorbed  to 
suspended  particulates,  and  dissolved  chemical  (McCay  and  Isaji,  2002).  The 
CHEMMAP  model  is  used  in  this  chapter  to  predict  the  propagation  of  chemicals. 

CHEMMAP  can  be  either  run  as  a  certain  scenario  with  specific  tidal  forcing, 
wind  and  spill  site  and  in  stochastic  mode  to  estimate  the  probable  distribution  and 
concentrations  resulting  from  hypothetical  spills.  The  next  chapter  used  six  different 
scenarios  including  two  locations  each  to  cover  all  possible  threats  and  different 
chemicals  in  San  Diego  Bay. 

The  model,  which  is  very  similar  but  not  the  same  as  WQMAP,  incorporates  a 
number  of  model  components  including  simulation  of  the  initial  release  for  surface  and 
subsurface  spills,  slick  spreading,  transport  of  floating,  dissolved  and  particulate 
materials,  evaporation  and  volatilization,  dissolution  and  adsorption,  sedimentation  and 
degradation. 

It  uses  physical  and  chemical  properties  such  as  density,  viscosity,  vapor  pressure, 
surface  tension,  water  solubility,  environmental  degradation  rates,  and 
adsorbed/dissolved  partitioning  coefficients.  CHEMMAP  relies  on  Stake’s  Law  in  order 
to  calculate  vertical  velocities.  Furthermore,  its  approach  towards  propagation  is 
Langrangian. 

The  outputs  of  the  model  include  the  trajectories,  and  concentrations.  More 
specifically,  it  is  possible  to  see  the  swept  area  by  a  floating  chemical,  as  well  as  the  total, 
absorbed,  dissolved  and  particulate  concentration  in  both  the  water  column  and  the 
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sediments.  The  most  important  is  that  it  is  then  possible  to  determine  the  range  of 
distances  and  directions  of  the  contamination  caused  from  the  spill  at  a  particular 
location. 

The  choice  of  the  chemicals  used  in  this  thesis  will  be  based  on  chemical/physical 
properties  and  toxicity  data  and  are  contained  in  a  database  compiled  from  published 
literature  sources,  mainly  French  et  al.  (1996)  and  Mackay  et  ah,  (1992  a,b,c,d).  Since 
several  properties  vary  with  temperature,  the  chemical  data  are  for  an  initial  temperature 
of  25°C.  The  model  corrects  these  parameters  to  the  ambient  temperature  for  the  spill 
incident.  The  algorithms  for  changing  viscosity  and  vapor  pressure  to  ambient 
temperature  are  taken  from  French  et  al.  (1996),  who  developed  regression  using  the  data 
in  Gambill  (1959).  For  pure  chemical  processes,  the  increase  per  10°C  is  assumed  to  be  a 
factor  of  2.  For  biological  processes  (e.g.,  degradation  rates),  the  increase  in  rate  per 
increase  of  10°C  is  assumed  to  be  a  factor  of  3  (McCay  and  Isaji,  2002). 

The  model  is  initialized  for  the  spilled  mass  at  the  location  and  depth  of  the 
release.  The  state  and  solubility  are  the  primary  determining  factors  for  the  initialization 
algorithm.  If  the  chemical  is  highly  soluble  in  water  and  is  either  a  pure  chemical  (e.g., 
the  benzene  scenario)  or  dissolved  in  water  (e.g.,  the  methanol  scenario),  the  chemical 
mass  is  initialized  in  the  water  column  in  the  dissolved  state  and  in  a  user-defined  initial 
volume.  For  insoluble  or  semi-soluble  gases  released  underwater  (e.g.,  the  naphthalene 
gas  scenario),  the  spilled  mass  is  initialized  in  the  water  column  at  the  release  depth  in  a 
user-defined  plume  volume,  as  bubbles.  The  median  particle  size  is  characterized  by  a 
user-defined  diameter  (McCay  and  Isaji,  2002). 

For  the  state  where  the  chemical  of  interest  is  both  adsorbed  to  particles  and 
dissolved  in  the  water  phase  of  the  bulk  liquid  (e.g.  our  ammonia  liquefied  gas  scenario), 
dissolved  mass  is  also  initialized  in  the  initial  plume  volume.  The  mass  of  chemical 
spilled  is  corrected  from  the  bulk  spill  volume  using  the  appropriate  density  and 
concentration  data  from  the  database  (McCay  and  Isaji,  2002). 

Chemical  mass  is  transported  in  three-dimensional  space  and  time,  by  surface 
wind  drift,  other  currents,  and  vertical  movement  in  accordance  with  buoyancy  and 
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dispersion.  The  model  simulates  adsorption  onto  suspended  sediment,  resulting  in 
sedimentation  of  material.  Stake’s  Law  is  used  to  compute  the  vertical  velocity  of  pure 
chemical  particles  or  suspended  sediment  with  adsorbed  chemical.  If  rise  or  settling 
velocity  overcomes  turbulent  mixing,  the  particles  are  assumed  to  float  or  settle  to  the 
bottom.  Settled  particles  may  later  re-suspend  (assumed  to  occur  above  20  cm/s  current 
speed).  Wind-driven  current  (drift)  in  the  surface  water  layer  (down  to  5m)  is  calculated 
within  the  fates  model,  based  on  hourly  wind  speed  and  direction  data.  Surface  wind  drift 
of  oil  has  been  observed  in  the  field  to  be  1-6%  of  wind  speed  in  the  direction  of  0-30 
degrees  to  the  right  (in  the  northern  hemisphere)  of  the  down- wind  direction  (Youssef 
and  Spaulding,  1993).  The  user  may  also  specify  the  wind  drift  speed  and  angle  (McCay 
and  Isaji,  2002). 

CHEMMAP  simulates  degradation,  volatilization,  evaporation,  dissolution, 
entrainment  and  spreading.  More  specifically,  spreading  is  simulated  using  the  Fay 
algorithm  (Fay,  1971).  Entrainment  is  modeled  as  for  oil  (Delvigne  and  Sweeney,  1988). 
Surface  floating  chemicals  interaction  with  shorelines  is  simulated  based  on  the 
algorithms  developed  for  oil  spills  (French  et  ah,  1999).  The  dissolution  rate  of  pure 
chemicals  is  a  function  of  solubility  using  a  first  order  constant  rate  equation.  The 
dissolved  chemical  in  the  water  column  is  assumed  to  adsorb  to  particulate  matter  in 
accordance  with  the  equilibrium  partitioning  theory  (DiToro  et  ah,  1991).  Evaporation  is 
modeled  following  the  theory  that  the  rate  of  mass  flux  to  the  atmosphere  increases  with 
vapor  pressure,  temperature,  wind  speed  and  surface  area  (Mackay  and  Matsugu,  1973). 
Volatilization  from  the  water  column  is  calculated  from  the  chemical’s  vapor  pressure 
and  solubility  (Fyman  et  ah,  1982).  Degradation  is  estimated  assuming  a  constant  rate  of 
“decay”  specific  to  the  environment  where  the  mass  exists  (i.e.,  atmosphere,  water 
column  or  sediment). 

The  spilled  chemical  is  modeled  using  the  Fagrangian  approach.  At  each  time 
step,  phase  transfer  rates  (evaporation,  dissolution,  volatilization,  and  entrainment)  are 
calculated  and  a  proportionate  percentage  of  the  spillets  are  transferred  to  the  new  phase 
(McCay  and  Isaji,  2002). 
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B.  ENVIRONMENTAL  DATA 

The  geographical  database  and  the  grid  used  for  WQMAP  are  also  used  in 
CHEMMAP  applications.  CHEMMAP  has  the  same  tidal  forcing  as  for  WQMAP,  and 
hourly  wind  data.  The  wind  is  more  influential  to  the  fate  of  insoluble  floating  chemicals 
than  contaminants  in  the  water  column. 

C.  CHEMICAL  ELEMENTS  DESCRIPTION 

In  the  current  threat  environment,  a  chemical  attack  in  a  big  city  hosting  a  large 
portion  of  the  U.S.  Naval  bases  (such  as  San  Diego)  is  anything  but  impossible.  Since  the 
thesis  is  unclassified,  it  is  not  permitted  to  use  elements,  which  would  be  used  in  the  case 
of  an  actual  weapon  with  mass  destruction  (WMD)  attack;  however  it  is  possible  to 
simulate  them  and  see  the  results  by  using  other  mainly  high  toxic  elements.  For  this 
reason,  the  choices  are  three  floating  chemicals  (methanol,  benzene  and  ammonia),  two 
sinking  ones  (chlorobenzene  and  trichloroethylene)  and  one  dispersing  in  the  air 
(naphthalene). 

1.  Floating  Chemicals 

Methanol8  (CH30H)  is  a  colorless  fairly  volatile  liquid,  belonging  to  the  aliphatic 
alcohols  chemical  type.  It  is  originally  distilled  from  wood,  but  currently,  it  is 
synthetically  produced  from  carbon  oxides  and  hydrogen.  It  has  a  faintly  sweet  pungent 
odor  like  that  of  ethyl  alcohol.  It  has  a  flash  point  12.222°C  and  its  density  is  791  kg/m3 
(at  25°C).  Its  vapors  are  slightly  heavier  than  air  and  may  travel  some  distance  to  a  source 
of  ignition  and  flash  back.  Any  accumulation  of  vapors  in  confined  spaces,  such  as 
buildings  or  sewers,  may  explode  if  ignited,  so  it  is  very  dangerous.  It  is  used  to  make 
chemicals,  to  remove  water  from  automotive  and  aviation  fuels,  as  a  solvent  for  paints 
and  plastics,  as  an  alternative  motor  fuel  and  as  an  ingredient  in  a  wide  variety  of 
products,  so  it  is  not  easy  control.  The  most  recent  inventory  is  estimated  to  be  1,125 
metric  tons  in  the  U.S.  alone. 

Methanol  reacts  violently  with  acetyl  bromide.  Mixtures  with  concentrated 
sulfuric  acid  and  concentrated  hydrogen  peroxide  can  cause  explosions.  It  reacts  with 
hypochlorous  acid  either  in  water  solution  or  mixed  water/carbon  tetrachloride  solution  to 

8  Other  synonyms  include  methyl  alcohol,  carbinol,  wood  alcohol,  wood  naptha  and  wood  spirit. 
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give  methyl  hypochlorite,  which  decomposes  in  the  cold  and  may  explode  on  exposure  to 
sunlight  or  heat.  It  gives  the  same  product  with  chlorine  and  can  react  explosively  with 
isocyanates  under  basic  conditions.  The  presence  of  an  inert  solvent  mitigates  this 
reaction.  A  violent  exothermic  reaction  occurred  between  methyl  alcohol  and  bromine  in 
a  mixing  cylinder.  A  flask  of  anhydrous  lead  perchlorate  dissolved  in  methanol  exploded 
when  it  was  disturbed. 

Hence,  methanol  is  high  flammable,  is  a  floater,  is  highly  volatile,  is  highly 
soluble  and  remains  dissolved.  Furthermore,  is  a  toxic  element,  and  therefore,  very 
dangerous.  Its  immediately  dangerous  to  life  or  health  indicator  (IDLH)  is  very  high 
(6,000  ppm),  its  short  tenn  exposure  limit  (STEL)  is  250  ppm  and  its  odor  threshold  is 
4.2  ppm.  Its  degradation  rate  is  0.097835  (%/day)  in  air  and  sediments  and  0.3024 
(%/day)  in  water. 

Even  though  it  is  not  “acutely”  toxic,  inhalation  can  cause  cough,  dizziness, 
headache,  nausea,  weakness  and  visual  disturbance.  Ingestion  can  be  even  worse  causing 
abdominal  pain,  convulsions,  shortness  of  breath,  unconsciousness  and  vomiting. 
However,  what  is  more  important  is  its  ecotoxicity.  For  certain  fish  (such  as  red  drums), 
shrimps,  mussels  and  snails,  it  can  be  lethal. 

Benzene9  (C6H6)  is  the  second  chemical  to  be  used.  It  is  also  toxic,  causing  the 
same  symptoms  as  methanol,  and  can  be  fatal  to  many  species  in  the  ecosystem  including 
all  the  aforementioned  ones  as  well  as  oysters,  clams,  trout,  salmon,  catfish  and  goldfish. 
As  regards  its  toxicity,  benzene  is  a  confirmed  carcinogen,  develops  and  reproduces 
toxins,  and  therefore,  is  extremely  dangerous.  The  water  maximum  contaminant  level  is 
0.000005  kg/m  .  Its  immediately  dangerous  to  life  or  health  indicator  (IDLH)  is  not  as 
high  as  methanol  (500  ppm),  its  short  term  exposure  limit  (STEL)  is  2.5  ppm  and  its  odor 
threshold  is  34  ppm  for  detection  and  97  ppm  for  recognition.  Its  degradation  rate  is 
0.97835  (%/day)  in  air,  0.097835  (%/day)  in  water  and  0.0097835  (%/day)  in  sediments. 

For  chemical  characteristics,  benzene  is  high  flammable,  floater,  highly  volatile, 
highly  soluble,  and  moderately  absorbable  to  particles.  It  is  a  clear  colorless  liquid  with  a 

9  Other  synonyms  are  benzol  and  cyclohexatriene. 
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petroleum-like  odor.  It  is  less  dense  than  water  (0.877  g/cm  )  and  slightly  soluble  in 
water.  Its  vapors  are  heavier  than  air.  Benzene  reacts  vigorously  with  alkyl  chloride  or 
other  alkyl  halides  even  at  -70°C  in  the  presence  of  ethyl  aluminum  dichloride  or  ethyl 
aluminum  sesquichloride  and  explosions  have  been  reported.  It  ignites  in  contact  with 
powdered  chromic  anhydride.  It  is  incompatible  with  oxidizing  agents  such  as  nitric  acid. 
Mixtures  with  bromine  trifluoride,  bromine  pentafluoride,  iodine  pentafluoride,  iodine 
heptafluoride  and  other  interhalogens  can  ignite  upon  heating. 

The  last  floating  chemical  used  is  ammonia10  (NH3)  liquefied  gas.  It  is  a  base.  It 
is  clear  and  colorless  and  has  a  strong  odor.  It  is  shipped  as  a  liquid  under  its  own  vapor 
pressure.  Its  density  in  the  liquid  fonn  is  12.8825  kg/m  .  Contact  with  the  unconfined 
liquid  can  cause  frostbite.  Gas  generally  regarded  as  nonflammable  but  does  bum  within 
certain  vapor  concentration  limits  and  with  strong  ignition.  Fire  hazard  increases  in  the 
presence  of  oil  or  other  combustible  materials.  Although  gas  is  lighter  than  air,  vapors 
from  a  leak  initially  hug  the  ground.  Prolonged  exposure  of  containers  to  fire  or  heat 
may  cause  violent  rupturing  and  rocketing.  Long-tenn  inhalation  of  low  concentrations 
of  the  vapors  or  short-term  inhalation  of  high  concentrations  has  adverse  health  effects. 
It  is  used  as  a  fertilizer  and  refrigerant,  and  in  the  manufacture  of  other  chemicals. 

Ammonia  is  a  floater,  highly  volatile,  highly  soluble  and  slightly  absorbable  to 
particles,  and  reacts  exothermically  with  all  acids.  Violent  reactions  are  possible.  It  also 
readily  combines  with  silver  oxide  or  mercury  to  form  compounds  that  explode  on 
contact  with  halogens.  When  in  contact  with  chlorates,  it  forms  explosive  ammonium. 

As  for  toxicity,  its  immediately  dangerous  to  life  or  health  indicator  (IDLH)  is 
relatively  small  (300  ppm),  its  short  term  exposure  limit  (STEL)  is  35  ppm  and  its  odor 
threshold  is  0.019  ppm.  Its  degradation  rate  is  0. 1586  in  both  air  and  water.  Contact  with 

ammonia  could  cause  skin  and  eye  bums  and  inhalation  some  burning  sensation,  cough, 
shortness  of  breath  and  sore  throat.  For  the  ecosystem,  its  slight  toxicity  can  be  lethal  to 
shrimp,  prawns,  salmon,  trout  and  catfish. 

2.  Sinking  Chemicals 

10  Other  synonyms  are  nitro-sil,  spirit  of  Hartshorn,  and  vaporole. 
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The  first  sinking  chemical  is  chlorobenzene11  (C6H5C1).  It  is  a  colorless  to  clear, 
yellowish  liquid  with  a  sweet  almond-like  odor.  It  is  insoluble  in  water  and  a  little  denser 
than  water  (1,107  kg/m  ).  Its  vapors  are  heavier  than  air.  It  is  used  to  make  pesticides, 
dyes,  and  other  chemicals,  chlorobenzene  undergoes  a  sometimes  explosive  reaction  with 
powdered  sodium  or  phosphorus  trichloride  and  sodium.  It  may  react  violently  with 
dimethyl  sulfoxide.  It  reacts  vigorously  with  oxidizing  agents.  It  attacks  some  forms  of 
plastic,  rubber  and  coatings  and  it  forms  a  shock  sensitive  solvated  salt  with  silver 
perchlorate. 

Chlorobenzene  is,  therefore,  a  sinker,  semi-volatile,  soluble,  highly  flammable 
and  moderately  absorbable  to  particles.  For  toxicity,  its  immediately  dangerous  to  life  or 
health  indicator  (IDLH)  is  quite  big  (1,000  ppm),  and  its  odor  threshold  is  1.3  ppm.  Its 
degradation  rate  is  0.09784  (%/day)  in  air,  0.0097835  (%/day)  in  water  and  0.00098 
(%/day)  in  sediments.  Inhalation  of  chlorobenzene  can  cause  drowsiness,  headache, 
nausea  and  unconsciousness.  Ingestion  causes  abdominal  pain.  As  regards  eco-toxicity,  it 
can  be  lethal  to  prawns,  trout  and  goldfish. 

The  second  sinker  is  trichloroethylene12  (C2HC13).  It  is  a  toxic  sinker  with  similar 
results  as  chlorobenzene  when  in  contact  with  humans.  Furthermore,  it  is  a  proven 
carcinogen.  It  is  a  clear  colorless  volatile  liquid  having  a  chloroform-like  odor.  It  is 
denser  than  water,  slightly  soluble  in  water  and  is  non-combustible.  It  is  used  as  a 
solvent,  fumigant,  in  the  manufacture  of  other  chemicals,  and  for  many  other  uses.  It  has 
been  determined  experimentally  that  mixtures  of  finely  divided  barium  metal  and  a 
number  of  halogenated  hydrocarbons  possess  an  explosive  capability.  Specifically, 
impact  sensitivity  tests  have  shown  that  granular  barium  in  contact  with 
monofluorotrichloromethane,  trichlorotrifluoroethane,  carbon  tetrachloride, 
trichloroethylene,  or  tetrachloroethylene  can  detonate.  It  has  been  determined 
experimentally  that  a  mixture  of  beryllium  powder  with  carbon  tetrachloride  or  with 
trichloroethylene  will  flash  or  spark  on  heavy  impact.  A  mixture  of  powdered  magnesium 

11  Other  synonyms  are  monochlorobenzene,  benzene  chloride  and  phenyl  chloride. 

12  Other  synonyms  are  TCE,  acetylene  trichloride,  algylen,  blacosolv,  chlorylen,  dow-tri,  ethylene 
trichloride,  fleck-flip,  tri-clene  and  trichloroethene. 
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with  trichloroethylene  or  with  carbon  tetrachloride  will  also  flash  or  spark  under  heavy 
impact. 

Thus,  trichloroethylene  is  a  sinker,  semi-volatile,  highly  soluble  and  moderately 
absorbable  to  particles.  As  regards  toxicity,  its  immediately  dangerous  to  life  or  health 
indicator  (IDLH)  is  the  same  as  chlorobenzene  (1,000  ppm),  its  short  term  exposure  limit 
(STEL)  is  100  ppm  and  its  odor  threshold  is  82  ppm.  Its  degradation  rate  is  0.09784 
(%/day)  in  air,  0.03024  (%/day)  in  water  and  0.003024  (%/day)  in  sediments.  Its  dangers 
for  the  eco-system  include  death  to  toads,  trout  and  flagfish. 

3.  Gaseous  Chemical 

The  last  chemical  is  naphthalene13  (gas)  (C10H8).  It  is  a  dark  liquid  mixture,  with 
much  different  qualities  than  all  the  previous  chemicals.  It  is  insoluble  in  water  and 
denser  than  water.  A  part  disperses  in  the  atmosphere  and  another  sinks  in  water.  Contact 
with  naphthalene  may  cause  irritation  to  skin,  eyes,  and  mucous  membranes.  It  can  cause 
confusion,  headache,  sweating,  nausea,  vomiting  and  jaundice  when  inhaled.  It  is  toxic 
by  ingestion  and  can  cause  abdominal  pain,  convulsions,  diarrhea,  dizziness  and 
unconsciousness. 

Its  toxicity  is  moderate  and  its  immediately  dangerous  to  life  or  health  indicator 
(IDLH)  is  250  ppm,  its  short  term  exposure  limit  (STEL)  is  15  ppm  and  its  odor  threshold 
is  0.038  ppm.  It  is  a  known  carcinogen,  and  therefore,  dangerous  to  humans.  A  mixture 
containing  naphthalene  may  react  vigorously  with  strong  oxidizing  agents.  It  can  also 
react  exothermically  with  bases  and  with  diazo  compounds.  Naphthalene  reacts  violently 
with  chromic  anhydride. 

Hence,  naphthalene  disperses  in  the  atmosphere  (but  can  be  also  a  sinker),  is  not 
volatile,  is  semi-to  non-soluble,  moderately  absorbable  to  particles,  is  highly  flammable 
and  does  not  react  rapidly  with  either  water  or  air.  Its  degradation  rate  is  0.97835  (%/day) 
in  air,  0.09784  (%/day)  in  water  and  0.00302  (%/day)  in  sediments.  Its  dangers  for  the 
eco-system  include  the  death  of  toads,  crabs,  shrimp,  cod,  salmon,  trout  and  oysters, 
being  one  of  the  most  dangerous  enemies  of  natural  underwater  life. 

13  Other  synonyms  are  naphthene,  tar  camphor  and  moth  balls. 
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Methanol 

Benzene 

Ammonia 

Chloro¬ 

benzene 

TCE 

Naphthalene 

(gas) 

Floatation 

Floater 

Floater 

Floater 

Sinker 

Sinker 

Sinker/  Air 
dispersed 

Solubility 

High 

High 

High 

Normal 

High 

Semi 

Volatility 

High 

High 

High 

Semi 

Semi 

None 

Absorption 

Dissolves 

Moderate 

Slight 

Moderate 

Moderate 

Moderate 

Flammability 

High 

High 

High 

High 

Water/Air 

rapid 

interaction 

No 

No 

No 

No 

Table  13.  Comparison  of  chemicals  used  in  our  CHEMMAP  scenarios. 
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VII.  TWO  REGIMES  OF  CHEMICAL  DISPERSION 


A.  METHANOL 

1.  Pollutants  Released  at  North  San  Diego  Bay 

Suppose  a  small  boat  drops  one  barrel  of  methanol  in  less  than  12  minutes  on 
midnight  July  4,  1993  (Independence  Day)  at  (32°43’N,  11 7°  13. 05’  W),  which  is  located 
in  the  northern  part  of  San  Diego  Bay.  The  release  depth  is  1  m  and  the  initial  plum 
thickness  is  0.5  m.  In  order  to  make  the  conditions  more  difficult  for  the  propagation, 
assume  that  there  is  no  wind  at  all. 

The  model  results  show  the  following  features:  In  three  hours,  the  methanol  is  in 
San  Diego  port  (Figure  30)  and  in  10  hours  it  is  spread  all  over  the  North  San  Diego  Bay 
(Figure  31).  In  16  hours,  it  reaches  the  Naval  Station.  However,  the  south  part  of  the  Bay 
is  contaminated  much  later.  After  two  days,  there  are  no  pollutant  particles  south  of 
32°40’N  (Figure  32).  After  three  days  the  Naval  Station  is  heavily  impacted  but  after  nine 
days,  there  are  still  no  pollutant  particles  south  of  32°39’N.  The  methanol  reaches  the 
south  end  of  the  Bay  only  after  20  days,  but  its  concentration  in  the  water  column  can  be 
neglected.  Figure  33  shows  the  final  swept  area  after  32  days. 

In  such  a  case,  it  can  be  concluded  that  there  is  plenty  of  time  to  take  protective 
measures  for  the  southern  part  of  the  Bay  where  the  results  of  such  an  incident  would  be 
minimal.  However,  for  the  port  of  San  Diego,  the  effect  is  immediate  and  serious.  For  the 
Naval  Station,  the  reaction  time  (19  hours  to  three  days)  is  critical. 
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Figure  30.  Methanol  dropped  in  North  San  Diego  Bay,  dissolved  concentration  in 

San  Diego  port/city  after  3  hours 


Furthermore,  after  five  entire  days,  one  third  of  the  methanol  is  still  in  the  water 
column  (Figure  34).  Note  that  it  takes  almost  12  days  for  the  concentration  in  the  water 
column  to  reach  10%  and  15  days  for  the  decayed  methanol  to  reach  a  level  of  80%. 
Moreover,  the  end-state  is  the  contamination  not  only  of  the  San  Diego  Bay  but  also  a 
considerable  part  of  the  sea  outside  the  Bay.  The  scenario  is  repeated  by  increasing  the 
amount  of  methanol,  but  nothing  changes  fundamentally.  The  mass  balance  curves  and 
the  area  contaminated  remain  the  same. 
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Figure  3 1 .  Methanol  dissolved  concentration  affecting  the  entire  northern  San  Diego 

Bay  in  10  hours. 
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Figure  34.  Mass  balance  for  methanol  dropped  in  North  San  Diego  Bay 
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2.  Pollutants  Released  at  South  San  Diego  Bay 

Suppose  a  small  boat  drops  one  barrel  of  methanol  in  less  than  12  minutes  on 
midnight  July  4,  1993  (Independence  Day)  at  (32°39’N,  1 17°07.92’  W),  which  is  located 
in  the  southern  part  of  San  Diego  Bay.  The  release  depth  is  1  m  and  the  initial  plum 
thickness  is  0.5  m.  In  order  to  make  the  conditions  more  difficult  for  the  propagation, 
assume  that  there  is  no  wind  at  all. 

The  model  results  show  the  following  features.  First  of  all,  very  few  pollutants 
reach  32°41’N  parallel  and  in  no  case  does  methanol  reach  San  Diego  port  (Figure  35). 
Flowever,  in  13  hours,  it  has  reached  the  Naval  Station  (Figure  36)  with  a  heavy  impact  in 
less  than  two  days.  In  any  case,  since  the  northern  part,  which  is  contaminated,  is  the 
Naval  Station,  it  is  safe  to  conclude  that  the  northern  San  Diego  Bay  is  more  important  as 
a  potential  target.  It  is  crucial  for  protective  measures  to  highlight  this  fact  because  a 
chemical  attack  in  the  South  San  Diego  Bay  will  have  minimal  effects,  or  at  least  much 
less  considerable  than  an  attack  (or  accident)  in  the  north  part  of  the  bay.  Figure  37  shows 
a  similar  but  different  result  as  regards  the  mass  balance  curves. 


Figure  35.  Total  swept  area  after  32  days  for  methanol  dropped  in  South  San  Diego 

Bay. 
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Figure  36. 


Methanol  dissolved  concentration  affecting  the  San  Diego  Naval  Station 

after  13  hours. 
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Figure  37.  Mass  balance  for  methanol  dropped  in  South  San  Diego  Bay 
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Thus,  the  decayed  methanol  reaches  80%  in  only  nine  days,  mainly  due  to  the 
inert  nature  of  methanol  in  combination  with  the  shallow  bathymetry  of  the  southern  part 
of  the  Bay.  It  is  important  to  single  out  that  in  the  first  case  (methanol  spill  over  in  the 
north),  the  dissolved  concentration  disappears  after  only  15  days,  but  in  the  second  case 
(south),  it  needs  29  days. 

The  final  conclusion  is  that  the  ecological  catastrophe  that  can  be  caused  with  a 
relatively  big  amount  of  methanol  spill  over  is  very  considerable,  especially  if  the  spill 
over  is  in  the  north.  It  can  also  be  harmful  to  humans. 

B.  BENZENE 

1.  Pollutants  Released  at  North  San  Diego  Bay 

Suppose  that  accidentally  or  on  purpose  there  is  a  spill  over  of  10  tons  of  benzene, 
whish  lasts  five  hours  and  happens  at  midnight  on  July  4,  1993,  first  in  the  same  location 
as  in  the  previous  case  in  the  northern  part  of  San  Diego  Bay.  Once  more,  the  release 
depth  is  1  m  and  the  initial  plum  thickness  is  0.5  m.  In  order  to  make  the  conditions  more 
difficult  for  the  propagation,  assume  again  that  there  is  no  wind. 

The  model  results  are  similar  to  the  methanol  scenario,  at  least  as  regards  the 
propagation:  In  a  little  more  than  three  hours,  the  benzene  reaches  San  Diego  port.  In  less 
than  12  hours,  there  is  already  a  heavy  impact  outside  the  bay  (Figure  38).  Nevertheless, 
the  southern  part  of  the  bay  is  contaminated  much  later.  The  Naval  Station  is  not 
contaminated  for  the  first  five  days  (Figure  39)  and  in  no  time  is  there  a  heavy  impact.  By 
the  time  the  benzene  reaches  the  Naval  Station  (5  days),  there  is  no  dissolved 
concentration  left.  After  eight  days,  there  is  no  particle  south  of  (p=  32  39N  and  the 
benzene  reaches  the  south  end  of  the  Bay  only  after  20  days,  when  its  percentage  in  water 
column  is  less  than  10%.  Again,  it  is  possible  to  conclude  that  in  such  a  case,  there  is 
plenty  of  time  to  take  protective  measures  for  the  southern  part  of  the  bay  (Chula  Vista 
area),  where  the  results  of  such  an  incident  would  be  minimal  (Figure  40). 

Furthermore,  starting  almost  immediately  and  lasting  for  the  entire  period,  a  large 
amount  of  benzene  (40%)  evaporates.  After  five  days  only,  20%  of  the  benzene  is  in  the 
water  column  (Figure  41).  However,  for  the  port  of  San  Diego,  the  effect  is  immediate 
and  serious,  as  seen  in  Figure  42. 
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Figure  38.  Benzene  dissolved  concentration  out  of  San  Diego  Bay  in  12  hours. 


Figure  39.  Benzene  swept  area  after  five  days. 
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Figure  40.  Benzene  swept  area  after  32  days. 
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Figure  41 .  Mass  balance  for  benzene  dropped  in  North  San  Diego  Bay. 
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Figure  42.  Absorbed  benzene  concentration  after  five  hours  in  North  San  Diego  Bay. 


The  decayed  percentage  never  reaches  more  than  20%.  Moreover,  the  end-state  is 
the  contamination  not  only  of  the  San  Diego  Bay  but  also  a  considerable  part  of  the  sea 
outside  the  Bay.  The  scenario  is  repeated  by  increasing  the  amount  of  benzene  and  the 
previous  results  verified.  The  mass  balance  curves  and  the  area  contaminated  remain  the 
same. 


2.  Pollutants  Released  at  South  San  Diego  Bay 

Suppose  that  accidentally  or  on  purpose  there  is  a  spill  over  of  10  tons  of  benzene, 
whish  lasts  five  hours  and  happens  at  midnight  on  July  4,  1993,  in  the  southern  part  of 
San  Diego  Bay,  in  a  location  closer  to  San  Diego  than  in  the  respective  methanol  scenario 
( cp  =32  39.34  N  and  A, =117  08.00W).  The  results  are  different,  but  again,  similar  to  the 
methanol  case.  First  of  all,  practically  benzene  particles  do  not  reach  San  Diego. 
Flowever,  the  Naval  Station  is  reached  within  12  hours  (Figure  43)  and  completely 
contaminated  in  less  than  three  days.  It  is  important  for  protective  measures  to  highlight 
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this  fact  because  a  chemical  attack  in  the  southern  part  of  the  Bay  would  affect  the  Naval 
Station,  and  in  the  northern  part,  would  affect  the  city.  Nevertheless,  after  17  days,  there 
is  dissolved  benzene  at  the  city/port  of  San  Diego  (Figure  44).  Figure  45  shows  the  end- 
state  of  benzene  swept  area. 


Figure  43.  Benzene  dissolved  concentration  at  the  San  Diego  Naval  Station  after  12 

hours. 
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Figure  44.  Benzene  dissolved  concentration  after  17  days. 


Figure  45.  Benzene  swept  area  after  32  days. 
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As  regards  the  mass  balance  curves,  there  is  a  similar  but  different  result,  as  seen 
in  Figure  46.  Comparing  Figures  41  and  46  shows  that  there  is  much  more  benzene  in  the 
water  column  (more  than  30%  after  10  days),  the  decayed  benzene  after  a  month  exceeds 
40%  and  only  the  rate  of  evaporation  is  the  same. 

The  final  conclusion  is  that  the  ecological  catastrophe  that  can  be  caused  with  a 
relatively  large  amount  of  benzene  spill  over  is  very  considerable  and  it  can  be  harmful  to 
humans  and  the  eco-system. 


Figure  46.  Mass  balance  for  benzene  dropped  in  South  San  Diego  Bay. 


C.  AMMONIA 

1.  Pollutants  Released  at  North  San  Diego  Bay 

Suppose  the  sinking  of  a  ship  containing  liquefied  ammonia  in  the  same  two 
locations  and  time  used  in  the  previous  scenario.  This  time  in  a  depth  of  3m  there  is 
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leakage  of  200  tons  of  liquefied  ammonia  which  is  completed  in  20  hours.  This  is  more 
of  an  environmental  threat  than  an  immediate  lethal  threat  to  humans.  Once  again,  there 
is  no  wind  forcing. 

The  results  are  similar  to  both  aforementioned  scenarios  as  regards  the 
propagation:  In  a  little  more  than  three  hours,  the  ammonia  reaches  San  Diego  port 
(Figure  47)  and  in  less  than  12  hours,  it  is  widespread  outside  of  the  bay  (Figure  48). 
Nevertheless,  the  southern  part  of  the  Bay  is  contaminated  much  later.  After  five  days, 
there  is  no  particle  south  of  cp  =  32  40N  and  the  benzene  reaches  the  southern  end  of  the 
Bay  only  after  15  days  (Figure  49).  The  Naval  Station  is  contaminated  in  a  period 
between  30  hours  (Figure  50)  and  five  days,  in  which  the  percentage  of  ammonia  in  the 
water  column  in  between  45  and  80%  (Figure  51).  After  15  days,  the  decayed  ammonia 
exceeds  50%  and  that  left  in  the  water  column  is  a  little  more  than  15%. 


It  is  again  possible  to  conclude  that  in  such  a  case,  there  is  plenty  of  time  to  take 
protective  measures  only  for  the  southern  part  of  the  Bay  (excluding  the  Naval  Station), 
where  the  results  of  such  an  incident  would  be  minimal.  Flowever,  for  both  the  city/port 
of  San  Diego  and  the  Naval  Station,  the  effect  is  immediate  and  serious.  The  end-state  is 
the  contamination  of  the  entire  San  Diego  Bay  (Figure  52)  as  well  as  a  considerable  part 
of  the  sea  outside  the  bay. 


Figure  47.  Ammonia  mean  dissolved  concentration  after  three  hours. 
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Figure  49.  Ammonia  mean  dissolved  concentration  after  15  days. 
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Figure  50.  Ammonia  mean  dissolved  concentration  after  30  hours. 
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Figure  5 1 .  Mass  balance  for  ammonia  dropped  in  North  San  Diego  Bay 
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Figure  52.  End-state  as  regards  swept  area  for  ammonia  dropped  in  North  San  Diego 

Bay 

2.  Pollutants  Released  at  South  San  Diego  Bay 

Suppose  the  sinking  of  a  ship  containing  liquefied  ammonia  in  the  southern  part 
of  the  Bay,  in  the  same  location  used  in  the  methanol  scenario  (cp=32  39.00N  and  X=1 17 
07.92W).  The  results  are  different  but  again  similar  to  the  previous  cases.  First  of  all, 
practically  benzene  does  not  damage  the  port  and  city  of  San  Diego  (Figure  53). 
Flowever,  in  12  hours,  it  has  reached  the  Naval  Station  (Figure  54)  and  in  less  than  three 
days,  it  has  a  very  heavy  impact  on  it.  It  is  important  for  protective  measures  to  highlight 
this  fact  because  a  chemical  attack  in  the  southern  part  of  the  bay  would  have  minimal 
effects  on  the  city,  which  will  be  affected  after  15  days  (Figure  55)  but  are  considerable 
at  the  Naval  Station. 

As  regards  the  mass  balance  curves,  there  are  similar  but  different  results,  as  seen 
in  Figure  56.  By  comparing  Figures  38  and  39,  notice  that  there  is  much  more  ammonia 
in  the  water  column  (more  than  50%  after  seven  days)  and  the  decayed  percentage  in  25 
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days  reaches  90%.  The  final  conclusion  is  that  the  ecological  catastrophe  that  can  be 
caused  with  a  relatively  large  amount  of  ammonia  spill  over  is  smaller  than  in  the 
previous  cases,  yet  considerable,  for  humans  and  the  eco-system. 


Figure  54.  Ammonia  mean  dissolved  concentration  after  12  hours. 
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Figure  55.  Ammonia  mean  dissolved  concentration  after  15  days. 
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Figure  56.  Mass  balance  for  ammonia  dropped  in  South  San  Diego  Bay. 
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D.  CHLOROBENZENE 


1.  Pollutants  Released  at  North  San  Diego  Bay 

Suppose  the  sinking  of  a  ship  containing  chlorobenzene,  in  the  same  two  locations 
and  time  used  in  our  previous  scenarios.  Again,  in  a  depth  of  3m,  there  is  leakage  of  200 
tons,  which  is  completed  in  20  hours.  Still,  there  is  no  wind  forcing.  The  results  are 
similar  to  all  aforementioned  scenarios  as  regards  the  propagation.  In  three  hours, 
chlorobenzene  reaches  San  Diego  port  (Figure  57)  and  in  less  than  12  hours,  it  is 
widespread  outside  of  the  bay  (Figure  58).  After  20  hours,  there  is  no  particulate 
concentration  and  after  30  hours,  the  Naval  Station  is  contaminated  (Figure  59). 
Nevertheless,  the  southern  part  of  the  bay  is  contaminated  much  later,  in  approximately 
19  days  (Figure  60),  and  finally  note  widespread  chlorobenzene  all  over  the  bay  (Figure 
61). 

Furthermore,  after  42  hours,  there  is  no  sediment.  There  is  very  little  decay  and  in 
no  case  exceeds  10%.  What  is  more  important  is  that  chlorobenzene  in  the  water  column 
decreases  from  almost  90%  to  45%  after  32  days,  which  is  a  very  impressive  percentage 
of  remaining  chlorobenzene  (Figure  62).  The  end-state  is  the  contamination  not  only  of 
the  San  Diego  Bay  but  also  a  considerable  part  of  the  sea  outside  the  bay. 


Figure  58.  Chlorobenzene  mean  dissolved  concentration  after  12  hours 


Figure  59.  Chlorobenzene  mean  dissolved  concentration  after  30  hours 
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Figure  60.  Chlorobenzene  mean  dissolved  concentration  after  19  days 
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2.  Pollutants  Released  at  South  San  Diego  Bay 

Suppose  the  sinking  of  a  ship  containing  chlorobenzene  in  the  southern  part  of  the 
Bay,  in  the  same  location  used  in  the  previous  scenarios  ( cp  =32  39.00N  and  A, =117 
07.92W),  the  results  are  very  similar  to  the  previous  cases.  First  of  all,  practically 
benzene  does  not  reach  the  port  and  city  of  San  Diego  (Figure  63).  However,  in  1 1  hours, 
it  has  reached  the  Naval  Station  (Figure  64). 

As  regards  the  mass  balance  curves,  there  is  a  similar  but  different  result,  as  seen 
in  Figure  65.  Comparing  Figures  62  and  65  shows  that  the  results  are  the  same 
concerning  sediments  and  particulate  concentrations.  However,  there  is  much  more 
chlorobenzene  in  the  water  column  (more  than  80%  after  32  days)  and  the  decayed 
percentage  reaches  15%.  The  final  conclusion  is  that  the  sinker  can  cause  a  greater 
ecological  catastrophe  that  a  floater  and  in  a  carefully  planned  attack  in  the  North  San 
Diego  Bay,  the  damage  might  be  overwhelming. 
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Figure  63.  Swept  area  of  chlorobenzene  dropped  in  South  San  Diego  Bay. 


Figure  64.  Chlorobenzene  reaches  the  Naval  Station  after  1 1  hours  when  dropped  in 

South  San  Diego  Bay 
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Figure  65.  Mass  balance  for  chlorobenzene  dropped  in  South  San  Diego  Bay. 


E.  TRICHLOROETHYLENE 

1.  Pollutants  Released  at  North  San  Diego  Bay 

Suppose  the  sinking  of  a  ship  containing  trichloroethylene,  in  the  same  two 
locations  and  time  used  in  the  previous  scenarios.  Again,  in  a  depth  of  3m,  there  is 
leakage  of  200  tons,  which  is  completed  in  20  hours.  As  in  all  previous  scenarios,  there  is 
no  wind  forcing.  The  results  are  not  surprisingly  the  same  as  regards  the  propagation.  In 
three  hours,  trichloroethylene  reaches  San  Diego  port  (Figure  66)  and  in  12  hours,  it  is 
widespread  outside  of  the  Bay.  A  small  difference  is  that  there  is  no  particulate 
concentration  and  a  very  small  percentage  of  sediments  even  after  32  days.  The  Naval 
Station  is  contaminated  after  30  hours  and  the  southern  part  of  the  Bay  is  contaminated  in 
approximately  19  days.  The  swept  area  end-state  is  very  similar  to  the  respective  case  of 
chlorobenzene  (Figure  67).  Furthennore,  the  decay  curve  reaches  25%  after  32  days. 
Trichloroethylene  in  the  water  column  remains  at  a  high  level  but  not  as  high  as 
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chlorobenzene  in  the  previous  scenario.  It  decreases  from  almost  90%  to  less  than  30%. 
Nevertheless,  the  first  12  days,  the  percentage  of  the  trichloroethylene  in  the  water 
column  is  above  50%  (Figure  68). 


Figure  67.  Swept  area  of  trichloroethylene  dropped  in  North  San  Diego  Bay. 
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Figure  68.  Mass  balance  for  trichloroethylene  dropped  in  North  San  Diego  Bay. 


2.  Pollutants  Released  at  South  San  Diego  Bay 

Suppose  the  sinking  of  a  ship  containing  trichloroethylene  in  the  southern  part  of 
the  Bay,  in  the  same  location  used  in  the  previous  scenarios  ( cp  =32  39.00N  and  A, =117 
07.92W),  the  results  are  very  similar  to  the  previous  cases.  Again,  trichloroethylene  does 
not  reach  the  port  and  city  of  San  Diego  (Figure  69).  However,  in  1 1  hours,  it  has  reached 
the  Naval  Station  (Figure  70). 

As  regards  the  mass  balance  curves,  there  is  a  similar  but  different  result,  as  seen 
in  Figure  71.  Comparing  Figures  68  and  71  shows  that  the  results  are  the  same 
concerning  sediments  and  particulate  concentrations.  However,  there  is  much  more 
trichloroethylene  in  the  water  column  (almost  60%  after  32  days)  and  the  decayed 
percentage  reaches  40%.  The  final  conclusion  is  that  the  sinker  can  cause  a  greater 
ecological  catastrophe  than  a  floater  but  is  much  less  significant  in  South  San  Diego  Bay 
than  in  the  northern  part  of  the  Bay. 
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Figure  69.  Swept  area  of  trichloroethylene  dropped  in  South  San  Diego  Bay. 


Figure  70.  Trichloroethylene  reaches  the  Naval  Station  after  1 1  hours  when  dropped 

in  South  San  Diego  Bay. 
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Figure  7 1 .  Mass  balance  for  trichloroethylene  dropped  in  South  San  Diego  Bay. 


F.  NAPHTALENE 

1.  Pollutants  Released  at  North  San  Diego  Bay 

Suppose  the  sinking  of  a  ship  containing  naphthalene,  in  the  same  two  locations 
and  time  used  in  the  previous  scenarios.  Again,  in  a  depth  of  3m,  there  is  leakage  of  200 
tons,  which  is  completed  in  10  hours.  Still,  there  is  no  wind  forcing.  Despite  the  much 
different  qualities  of  this  chemical,  the  results  are  identical  to  the  aforementioned 
scenarios  as  regards  the  propagation.  In  three  hours,  naphthalene  gas  reaches  San  Diego 
port  (Figure  72)  and  in  less  than  12  hours,  it  is  widespread  outside  of  the  Bay.  After  30 
hours,  the  Naval  Station  is  contaminated.  The  southern  part  of  the  Bay  is  contaminated 
much  later,  in  approximately  20  days  and  finally,  widespread  naphthalene  is  observed  all 
over  the  bay  (Figure  73).  As  regards  the  mass  balance  curves,  the  final  decay  is  a  bit 
more  than  30%.  After  the  fourth  day,  the  naphthalene  in  the  water  column  is  around  35%, 
and  after  32  days,  it  is  only  5%.  After  1 1  days,  both  the  decay  and  the  water  column 
percentage  are  at  23%  (Figure  74).  The  most  interesting  difference  is  the  existence  of 
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important  absorbed  concentration.  After  18  hours,  there  is  much  concentration  in  the 
port/city  of  San  Diego  (Figure  75)  and  after  nine  days,  there  is  still  some  at  the  mouth  of 
the  bay  (Figure  76). 


Figure  72.  Naphthalene  mean  dissolved  concentration  after  three  hours. 
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Figure  73.  Swept  area  after  32  hours  of  naphthalene  dropped  at  North  San  Diego 

Bay. 
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Figure  74.  Mass  balance  for  naphthalene  dropped  in  North  San  Diego  Bay. 


Figure  75.  Naphthalene  mean  absorbed  concentration  after  18  hours. 
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Figure  76.  Naphthalene  mean  absorbed  concentration  after  nine  days. 

2.  Pollutants  Released  at  South  San  Diego  Bay 

Suppose  the  sinking  of  a  ship  containing  naphthalene  in  the  southern  part  of  the 
Bay,  in  the  same  location  used  in  the  previous  scenarios  ( cp  =32  39.00N  and  A, =117 
07.92W),  the  results  are  very  close  to  the  previous  cases.  Naphthalene  does  not  reach  the 
port  and  city  of  San  Diego  (Figure  77).  Flowever,  in  12  hours,  it  has  reached  the  Naval 
Station,  where  there  is  also  an  important  absorbed  concentration  (Figures  78  and  79). 
Moreover,  absorbed  concentration  is  still  high  after  five  days  and  six  hours  (Figure  80). 

As  regards  the  mass  balance  curves,  there  is  a  different  result,  as  seen  in  Figure 
81.  The  decay  curve  is  exactly  opposite  to  the  water  column  one,  and  after  12  days,  they 
are  both  at  50%.  The  end-state  for  the  water  column  is  a  little  less  than  20%  and  the 
decayed  naphthalene  exceeds  80%.  The  final  conclusion  is  that  the  semi  to  non-soluble 
naphthalene  can  cause  a  very  large  ecological  catastrophe,  probably  more  than  any  other 
sinker  and  definitely  more  than  the  floaters.  Still,  it  is  much  less  significant  in  South  San 
Diego  Bay  than  in  the  northern  part  of  the  bay. 
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Figure  77.  Swept  area  after  32  hours  of  naphthalene  dropped  at  South  San  Diego 

Bay. 


Figure  78.  Naphthalene  mean  dissolved  concentration  after  12  hours. 
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Figure  80.  Naphthalene  mean  absorbed  concentration  after  126  hours. 
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Mass  Balance  for  Naphthalene  (gas) 


T  — 

Evaporated 


Figure  8 1 .  Mass  balance  for  naphthalene  dropped  in  South  San  Diego  Bay. 


G.  OVERALL  ASSESSMENT 

An  overall  assessment  is  that  San  Diego  Bay,  being  a  semi-closed  tidal  basin,  is  a 
probable  target  for  a  terrorist  chemical  attack  or  a  location  where  a  chemical  accident 
could  cause  severe  problems.  The  sensitive  eco-system  would  be  impacted  both  inside 
and  outside  the  Bay.  If  the  chemical  is  a  sinker,  the  results  are  more  catastrophic  than  if  it 
were  a  floater.  Since  the  water  exchange  with  the  Pacific  Ocean  occurs  only  through  a 
narrow  entrance,  the  water  would  be  contaminated  for  long  time. 

The  propagation  model  shows  that  the  northern  part  of  the  bay  is  more  likely  to  be 
a  target  because  it  would  affect  the  city  and  the  Naval  Station,  and  it  would  reach,  even 
slightly,  the  South  San  Diego  Bay  and  would  spread  outside  the  bay  as  well.  Any  attack 
or  accident  in  the  southern  part  (south  of  the  Naval  Station)  would  have  minimal  effects 
on  the  city  and  the  shores  of  Coronado  Island  and  none  outside  the  bay.  On  the  other 
hand,  when  the  spill  occurs  in  the  southern  part  of  the  bay,  a  larger  percentage  of  the 
chemical  remains  in  the  water  column  and  for  longer  period  of  time. 
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In  3  hours:  San  Diego  port/city 

In  10  hours:  Entire  North  SD  Bay 

In  12  hours:  Outside  SD  Bay 

In  16-30  hours:  Naval  Station 

In  5  days:  Heavy  impact  on  Naval  Station 

In  20  Days:  South  end  of  the  Bay 

In  32  Days:  The  entire  SD  Bay 

Table  14.  Results  in  North  San  Diego  Bay 
In  12  hours:  Naval  Station 

In  15-17  days:  Small  part  of  absorbed  or  dissolved  chemical  in  San 
Diego  city/port 

After  32  days:  No  effect  to  North  San  Diego  Bay 

Table  15.  Results  in  South  San  Diego  Bay 


Methanol:  After  3  days,  45-50%  in  water  column,  after  20  days,  less  than  5% 
-  rest  decayed. 

Benzene:  45%  evaporates.  After  2  days,  30-50%  in  water  column,  after  20 
days  8-18%  -  rest  decayed. 

Ammonia:  After  3  days,  50-75%  in  water  column,  after  20  days,  8-18%  -  rest 
decayed 


Table  16.  Results  for  floaters  in  North  San  Diego  Bay  (first  percentage)  and  South 

San  Diego  Bay  (second  percentage). 
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Chlorobenzene:  After  5  days,  65  -  97%  in  water  column,  after  20  days,  50- 
90%  -  rest  decayed. 

Trichloroethylene:  After  5  days,  60-93%  in  water  column,  after  20  days,  38- 
71%  -  rest  decayed. 

Naphthalene  (gas/air  dispersed):  After  5  days,  33  -  78%  in  water  column,  after 
20  days,  12-33%  -  rest  decayed. 


Table  17.  Results  for  sinkers  in  North  San  Diego  Bay  (first  percentage)  and  South 

San  Diego  Bay  (second  percentage). 


Taking  into  consideration  that  the  database  available  does  not  contain  highly 
acute  toxic  elements  or  any  elements  that  consist  of  a  chemical  threat,  it  is  necessary  to 
highlight  the  vulnerability  of  the  Bay  in  case  of  such  an  attack  due  to  its  oceanography/ 
topography  and  bathymetry. 

Results  concerning  San  Diego  Bay  can  also  be  applied  to  studies  in  other  semi- 
closed,  barotropic,  no-wind  driven  circulation  basins. 
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VIII.  CONCLUSIONS  -  RECOMMENDATIONS 


This  study  determined  that  a  2D  depth-integrated  model  can  satisfy  certain  Navy 
applications  in  coastal  bays.  In  order  to  arrive  to  such  a  conclusion  the  thesis  first  showed 
that  San  Diego  Bay  is  barotropic.  Then,  as  regards  the  forcing  mechanisms,  it  proved  that 
the  wind  is  not  a  significant  factor.  Hence,  in  a  tidal  semi-enclosed  basin  like  San  Diego  a 
2-D  model  can  be  used  for  time  scales  ranging  from  “instantaneous”  (what  a  diver  will 
feel  in  the  water)  to  days  (the  drift  of  floating  objects).  However,  we  should  notice  that  in 
the  entrance  of  a  semi-enclosed  bay  we  should  be  more  careful  when  applying  a  2-D 
model  (depending  on  the  accuracy  we  wish  to  achieve),  because  both  the  influence  of  the 
wind  and  the  system  of  the  currents  outside  the  bay  cannot  be  ignored. 

As  regards  the  evaluation  of  the  two  models  purchased  by  U.S.  Naval 
Oceanographic  Office,  the  conclusion  is  that  they  need  improvement.  The  main 
consideration  for  applying  them  in  littoral  waters  should  be  the  smoothing  of  the 
bathymetry.  On  the  one  hand  it  is  essential  for  the  model  to  run  properly  and  on  the  other 
hand  might  produce  significant  errors  in  very  shallow  areas  (like  the  southern  part  of  San 
Diego  Bay). 

Finally,  the  study  shows  the  vulnerability  of  a  semi-enclosed  tidal  basin  in  a 
possible  chemical  attack  or  accident,  with  the  aforementioned  particular  results  for  San 
Diego  Bay.  In  order  to  summarize  these  results,  it  should  be  repeated  that  in  a  case  of  a 
chemical  attack  or  accident,  first  the  sensitive  eco-system  would  be  severely  damaged,  no 
matter  the  nature  of  the  event  and  the  location.  If  the  chemical  were  a  sinker,  the  results 
would  be  more  catastrophic  than  if  it  were  a  floater.  Since  the  water  exchange  with  the 
Pacific  Ocean  occurs  only  through  a  narrow  entrance,  the  water  would  be  contaminated 
for  long  time. 

Two  regimes  of  the  chemical  dispersion  were  found  in  this  thesis.  The  first  was 
the  case  of  an  attack/accident  in  the  North  San  Diego  Bay.  In  that  case  the  entire  bay 
would  be  contaminated.  In  3  hours  the  chemical  would  reach  San  Diego  port  and  city,  in 
12  hours  the  entire  northern  part  of  the  bay  would  be  affected  and  in  2-5  days  the  Naval 


91 


Station  (located  in  the  south  part  of  the  bay)  would  be  contaminated  as  well.  The  rest  of 
the  Bay  would  be  reached  much  later.  The  second  regime  was  an  attack/accident  in  the 
South  San  Diego  Bay.  In  such  case,  the  incident  would  have  minimal  effects  on  the  city 
and  the  shores  of  Coronado  Island  (located  in  the  north  part  of  the  bay)  and  none  outside 
the  Bay.  On  the  other  hand,  when  the  spill  occurs  in  the  southern  part  of  the  Bay,  a  larger 
percentage  of  the  chemical  remains  in  the  water  column  and  for  longer  period  of  time, 
which  makes  it  more  “effective”,  which  in  a  case  of  a  chemical  attack  means  lethal. 

For  the  aforementioned  reasons,  the  propagation  model  shows  that  the  northern 
part  of  the  bay  is  more  likely  to  be  a  target  because  it  would  affect  the  city  and  the  Naval 
Station,  and  it  would  reach,  even  slightly,  the  South  San  Diego  Bay  and  would  spread 
outside  the  bay  as  well.  In  general,  results  concerning  San  Diego  Bay  can  also  be  applied 
to  studies  in  other  semi-closed,  barotropic,  no-wind  driven  circulation  basins. 

As  regards  recommendations  for  future  research,  it  should  be  mentioned  that  the 
use  of  more  accurate  bathymetry  and  of  a  finer  grid  would  give  better  results  in  a  similar 
case.  Moreover,  the  use  of  more  recent  ADCP  measurements,  during  a  longer  period  of 
time  would  further  improve  the  results  and  verify  the  overall  conclusions.  It  would  be 
helpful  if  the  ADCPs  used  in  the  future  were  located  in  a  bigger  distance  from  the  shore. 

A  more  detailed  comparison  of  3D  vs.  2D  model  is  encouraged,  as  well  as  its 
application  for  drift  and  for  instantaneous  current  prediction.  Last  but  not  least,  as  regards 
chemical  propagation,  a  classified  research  with  data  unavailable  to  foreigners  about  real 
chemical  threats  (e.g.  anthrax)  should  be  conducted. 
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